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We review the currently available formalisms for radiative energy loss of a high-momentum parton 
in a dense strongly interacting medium. The underlying theoretical framework of the four commonly 
used formalisms is discussed and the differences and commonalities between the formalisms are 
highlighted. A quantitative comparison of the single gluon emission spectra as well as the energy 
loss distributions is given for a model system consisting of a uniform medium with a fixed length 
of L = 2 fm and L — 5 fm (the 'Brick'). Sizable quantitative differences are found. The largest 
differences can be attributed to specific approximations that are made in the calculation of the 
radiation spectrum. 



I. INTRODUCTION 

QCD jets are produced in hadronic interactions of all 
kinds, by the hard scattering of the constituent partons 
(quarks and gluons) of the colliding projectiles. One of 
the most significant experimental discoveries in the col- 
lision of heavy nuclei at -^/saT/v = 200 GeV at the Rela- 
tivistic Heavy Ion Collider (RHIC) was the observation 
of a strong suppression of the inclusive yield of high mo- 
mentum (high px) hadrons p] E] and the semi- inclusive 
rate of azimuthally back-to-back high px hadron pairs 
relative to expectations from p+p and d+Au collisions 
[5HS], as expected from jet quenching by the hot and 
dense medium formed in these collisions. Evidence for 
jet quenching has also been observed recently in Pb+Pb 
collisions at the LHC [SHE]- 

Most generally, the term 'jet quenching' refers to the 
modification of the evolution of an energetic parton in- 
duced by its interactions with a colored medium. The 
expected manifestations of modified evolution of parton 
showers include the suppression of hadronic spectra at 
high pt and their back-to-back azimuthal correlations, 
as well as the enhancement of hadron spectra at low px 
and possibly the angular broadening of internal jet struc- 
ture and di-jet acoplanarity. 

Data from the RHIC heavy ion program strongly sup- 
port the picture that jet quenching is caused by the loss 
of energy of the primary parton, either by collisions with 



constituents of the medium (collisional or elastic energy 
loss [HI HHj) or by gluon bremsstrahlung (radiative or 
inelastic energy loss |11H13| ). prior to hadronization in 
the vacuum. The loss of energy suffered by an energetic 
quark or gluon penetrating a QCD medium probes dy- 
namical properties of the medium. Jet quenching pro- 
vides unique and powerful tools for studying the proper- 
ties of the hot and dense matter produced in heavy ion 
collisions. 

In recent years, a number of different dynamical mod- 
els of jet quenching has been formulated and compared to 
the jet quenching signatures measured at RHIC (see [TH - 
I17j for reviews). These models are based on the common 
assumption that interactions of the energetic parton and 
the radiated gluon with the medium can be calculated via 
perturbative QCD (pQCD), regardless whether the prop- 
erties of the medium itself can be treated perturbatively. 
Within such perturbative approaches, however, different 
approximations exist in the calculation of radiative en- 
ergy loss of the leading parton. Direct comparisons of 
the various pQCD-based models show large quantitative 
discrepancies in the medium density that is needed to de- 
scribe the same suppression of the inclusive production 
of light leading mesons [T8H2U] . The full modeling of jet 
quenching in the complex, dynamic process of heavy ion 
collisions contains a large number of components, includ- 
ing the initial production spectrum, a time-dependent 
medium density profile, and fragmentation. The com- 
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plexity of the modeling, with each model at present em- 
ploying a different dynamical framework, makes it diffi- 
cult to isolate the specific differences between the energy 
loss formalisms, and thereby to constrain the underlying 
mechanisms of jet quenching through full model compar- 
isons. 

In this paper we review the current state of the art 
of pQCD-based radiative energy loss models, including 
a discussion of their theoretical foundations and limita- 
tions, their kinematic region of applicability, as well as 
various model-dependent sources of uncertainty. In or- 
der to make a direct comparison of the models, we pro- 
vide a systematic quantitative comparison of energy loss 
formalisms using a highly simplified model problem: the 
energy loss of light quarks in a static medium of constant 
density and a fixed length, known as the QGP brick. We 
will restrict the discussion to leading parton energy loss, 
for which the largest body of theoretical work exists. 

The paper is organized as follows. The first section pro- 
vides an executive summary of our main findings, with 
particular emphasis on common features, technical and 
conceptual differences, and uncertainties of different jet 
quenching models. This section is a self-contained nar- 
rative with a short introduction to the different classes 
of jet quenching models. It is intended to be accessi- 
ble to a non-expert reader. Many of the technical details 
will be omitted from the executive summary, which refers 
to following sections for more detailed explanations and 
numerical studies supporting the main statements. For 
more extensive descriptions of the various jet quenching 
formalisms and of the available phenomenology including 
model-data comparisons, we refer to the existing reviews 



II. EXECUTIVE SUMMARY 

A. Jet quenching is caused by partonic interactions 
in dense matter 

We first discuss the factorized perturbative QCD 
framework for calculating jet and hadron production in 
the absence of medium effects. In elementary collisions 
(e + e~, pp, and pp), our understanding of hadron pro- 
duction at high transverse momentum is relatively ma- 
ture. High-p^ hadron production results from partonic 
processes with large momentum transfer that can be de- 
scribed with controlled uncertainty as the convolution of 
incoming parton distribution functions (PDFs), a hard 
partonic collision process, and the fragmentation of the 
partonic final state [21] . The PDFs and Fragmentation 
Functions (FFs) are universal, non-perturbative, scale- 
dependent distributions that obey DGLAP evolution and 
are determined by global fits to data from elementary col- 
lisions 

In ultra-relativistic heavy ion collisions, there is cur- 
rently no firm theoretical argument that production cross 
sections factorize. Rather, factorization is a working as- 



sumption that is consistent with phenomenological anal- 
yses made so far and that underlies all models discussed 
in this report. General considerations indicate that if 
factorization is assumed, the hard partonic interaction 
itself cannot be modified by the medium, since it occurs 
on temporal and spatial scales too short to be resolved 
by the medium. However, it is also expected on gen- 
eral grounds that the parton distribution functions in 
nuclei differ from those in nucleons, and it is possible 
that the fragmentation of outgoing partons is modified 
by the presence of a hot QCD medium. 

We know with certainty that jet quenching receives its 
main contribution on the level of the outgoing fragment- 
ing parton. This follows from several lines of argument. 
In particular, it is found that the production of direct 
photons, which have no strong interaction with medium, 
does not show the large suppressions characteristic for 
jet quenching |27j . Moreover, there are by now several 
global PDF fits that parametrize the nuclear dependence 
of incoming parton distribution functions (nPDFs). For 
the typical momentum fractions x and virtualities Q rel- 
evant for hard processes, these nPDF-fits show nuclear 
modifications of the order of 10-30 % that cannot ac- 
count for the observed factor 4-5 suppression of single 
inclusive hadron spectra [28l [29] . Therefore, the energy 
degradation observed in nucleus-nucleus collision for all 
hadronic high transverse momentum spectra must be due 
to dynamical effects occurring after the hard process ('fi- 
nal state effects'). 

Furthermore, there is strong experimental evidence 
that the final state effect responsible for jet quenching is 
of partonic nature, i.e., that it occurs prior to hadroniza- 
tion. In particular, correlations of high px hadron pairs 
with small angular separation are largely unmodified in 
systems that exhibit strong inclusive suppression [3H5]- 
In addition, there is no evidence that the characteristi- 
cally different nuclear 'absorption' cross sections of dif- 
ferent hadron species play a role in understanding the jet 
quenching effect. 

As a consequence of these generic observations, all 
current efforts to understand and simulate jet quench- 
ing are based on models of parton energy loss. Fur- 
ther progress then depends on promoting the relation be- 
tween jet quenching phenomenology and parton energy 
loss from a qualitative to a more and more quantitative 
one. To this end, we embark in the following on a critical 
assessment of current parton energy loss models. 

B. Models of radiative parton energy loss 

The elementary branching process within a medium of 
finite size L is depicted in Fig. [T] Calculations of ra- 
diative parton energy loss aim at determining such pro- 
cesses within the framework of perturbative QCD. All 
approaches studied so far involve significant assumptions 
and approximations about 

• the virtuality and (repeated) branching of the hard 
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FIG. 1: Schematic view of a hard parton-parton interactions 
(grey blob) producing a highly energetic parton, which subse- 
quently undergoes parton branching processes. In a heavy ion 
collision, this parton evolution occurs within a dense medium 
(blue area) that interferes with the vacuum evolution via a 
priori unknown elastic and inelastic interactions. 



parton; 

• the nature of the medium through which the ener- 
getic parton propagates; and 

• kincmatical approximations for the interaction be- 
tween medium and projectile parton. 

It will turn out that the largest quantitative differences 
among the various models arise not so much from dif- 
ferences in their basic underlying assumptions about the 
medium and its interaction with the high energy par- 
ton, but rather in the implementation of the simplifying 
approximations made to carry out the derivations and 
calculations, most importantly, the treatment of energy- 
momentum constraints and large angle radiation. 

In the remainder of this Section, we will first provide 
a short general discussion of in-medium QCD radiation, 
before describing the specific implementations. 



1. Virtuality and parton branching in the medium 

A high-pr parton, produced by a hard initial collision 
between incoming partons, carries initially a high virtu- 
ality. Even in the absence of a medium, the parton will 
undergo 'vacuum' splitting processes to reduce its off- 
shellness. When the radiations occurs in the medium, the 
question arises whether this parton splitting is just the 
same as in the vacuum, or rather an additional medium- 
induced splitting, or a result of the interference between 
medium-induced effects and dynamics that would also 
occur in the absence of a medium. None of the existing 
parton energy loss calculations treats the entire dynamics 
depicted in Fig. [I] in a field-theoretically rigorous fashion. 

The medium effect on parton splitting is brought about 
by interactions of the high-energy parton and the radi- 
ated gluon with the medium. For small-angle scattering, 



the parton and the radiated gluon propagate along sim- 
ilar paths, leading to significant interference and a finite 
formation time of the gluon, which suppresses the gluon 
radiation compared to incoherent emission. This effect 
is often referred to as the Landau-Pomeranchuk-Migdal 
effect. A characteristic consequence of this interference is 
that the amount of energy loss grows quadratically with 
L for in-medium path lengths that are small compared 
to the formation time [T2] . 

In general, multiple splittings may occur in the 
medium. A full calculation of multi-gluon final states 
would include interference terms between the different 
emitted gluons. Recently some progress has been made 
on this particular problem, using calculations of the in- 
terference between two emitters (the antenna problem), 
which show interesting features as angular ordering or 
partial decoherence [30l433j . The generalization of these 
results to multiple gluon emission is not known at present 
and all current calculations are based on repeated ap- 
plication of a single-gluon emission calculation. In the 
following we will separately discuss the single-gluon radi- 
ation kernel and the prescription used for multiple gluon 
emission. 



2. Modeling the medium 

The main interest in studying jet quenching is to use it 
to characterize the medium through which the projectile 
parton propagates. It is of great interest to implement 
different models of the medium in jet quenching calcula- 
tions to determine how the medium model and properties 
affect jet quenching. In the current practice, each model 
description uses a particular set of simplifications, which 
can be classified as follows: 

1. The medium is modeled as a collection of static 
scattering centers 

In this approach, the medium is modeled as a set 
of static colored scattering centers with a specified 
density distribution along the trajectory of the pro- 
jectile. A decreasing density approximates the ef- 
fects of an expanding medium. Because gluon ra- 
diation from the scattering centers is ignored, cal- 
culations in this set-up lead to gauge-invariant (i.e. 
physically meaningful) results only to leading order 
in a high-energy approximation. By construction 
this set-up neglects recoil effects and thus does not 
allow for elastic parton energy loss. 

This medium model was pioneered by Baier, Dok- 
shitzer, Mueller, Peigne and Schiff [T21 IM1 [33] 
(BDMPS) and independently by Zakharov [36 ] [37 ] . 
Gluon radiation is formulated in a path-integral 
that resums scatterings on multiple static colored 
scattering centers. Wiedemann [38] showed how 
this path-integral can be used to include interfer- 
ence effects between vacuum and medium-induced 
radiation in such a way that also the fc^-differential 
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medium-induced gluon distribution is accounted 
for. 

In the original BDMPS derivation of the medium- 
induced gluon distribution [12j [34] , the soft gluon 
approximation, x <C 1, was used. In later deriva- 
tions [35], the radiation spectrum is multiplied by 
an overall splitting function to take corrections for 
finite x into account. In the literature, the term 
"BDMPS-Z formalism" has been used for both ver- 
sions of the formalism. The difference between the 
radiation spectrum with and without this splitting 
function P s _> g (x) is discussed in Section III C 5 and 
Fig. [14] 

Most analytical and numerical results from the 
BDMPS-Z formalism use a saddle point approxi- 
mation that amounts effectively to assuming that 
the projectile interacts with the medium via multi- 
ple soft scattering processes. In the following, the 
numerical results from the low-cc, multiple soft scat- 
tering implementation of the BDMPS-Z formal- 
ism are based on the work by Armesto, Salgado 
and Wiedemann [55H4T?] . abbreviated as ASW-MS. 
In the totally coherent limit, in which the entire 
medium acts coherently towards gluon production, 
the multiple-soft scattering formalism results in a 
radiation spectrum that is a radiation term for 
gluon production with momentum transfer con- 
voluted with a Gaussian elastic scattering cross 
section oc exp [— q^_/q L] . In this limit, the 
medium is fully characterized by the transport co- 
efficient <J, the mean of the squared transverse mo- 
mentum exchanged per unit path length. 

The opacity expansion was pioneered by Gyulassy, 
Levai and Vitev [4TJ [42] (GLV) and independently 
by Wiedemann [35]. It also includes the interfer- 
ence between vacuum and medium-induced radi- 
ation and is based on a systematic expansion of 
the calculation in terms of the number of scatter- 
ings. Since the BDMPS-Z path-integral formalism 
can serve as a generating functional for the opac- 
ity expansion [3S] , the opacity expansion formalism 
is another limit for the solution of the BDMPS-Z 
path integral. 

In most existing calculations, only the leading term 
(N = 1) is included, but the behaviour for larger 
opacities has been explored in [43 , 44J . The medium 
is characterized by two model parameters, the den- 
sity of scattering centers n or mean free path A, and 
a Debye screening mass hd used to regulate the 
infrared behavior of the single scattering cross sec- 
tion. In contrast to the multiple soft scattering ap- 
proximation, this approach includes the power-law 
tail of the scattering cross section expected from 
QCD, leading to shorter formation times of the ra- 
diation compared to the multiple-soft scattering ap- 
proximation |45| . There are several model imple- 
mentations in the literature that differ significantly 



in their kinematic approximations. We will use two 
implementations (ASW-SH and DGLV (46J [4 7]) to 
discuss the differences in detail in Section IIII Bl 
Extension of the opacity expansion formalism to 
a medium with dynamical scattering centers was 
recently explored by Djordjevic and Heinz [48] , 

2. The medium is characterized by matrix elements of 
gauge field operators 

In principle, multiple gluon exchanges between 
a partonic projectile and a spatially extended 
medium can be formulated in a field theoretically 
rigorous fashion by describing the medium in terms 
of expectation values of 2-, 4-, 6-, 8-, ... field corre- 
lation functions. Energy loss calculations based on 
the higher-twist (HT) approach were pioneered by 
Guo and Wang ISO] • The approach includes the 
interference between vacuum and medium-induced 
radiation. Properties of the medium enter the cal- 
culation in terms of higher-twist matrix elements. 
In practice the matrix elements are factorized in the 
nPDF and matrix elements describing the interac- 
tion between final state partons and the medium. 
This factorization is valid at leading order in the 
path length L in the medium. As we discuss in 
some detail in Section HID the approximations 



currently employed result in a formulation of par- 
ton energy loss calculations that closely resembles 
models starting from a set of static scattering cen- 
ters. 

3. Thermally equilibrated, perturbative medium 

The formulation in rigorous field-theory of parton 
energy loss in a weakly-coupled medium in per- 
fect thermal equilibrium was developed by Arnold, 
Moore and Yaffe [52 [5J] (AMY). The medium is 
formulated as a thermal equilibrium state in Hard 
Thermal Loop improved finite temperature pertur- 
bation theory. As a consequence, all properties 
of the medium are specified fully by its temper- 
ature and baryon chemical potential. The calcu- 
lation does not incorporate vacuum branching of 
the projectile parton. In principle, the perturbative 
description of the thermal medium applies only at 
very high temperature T ^> T c . 

As seen from the list above, different models of parton 
energy loss characterize the medium in terms of different 
primary model parameters. In the existing literature, 
the different approximations used for the medium in the 
various approaches have led to different ways to specify 
the medium properties. Recently it has become custom- 
ary to translate the primary model parameters into an 
effective q, that has the physical interpretation of an av- 
eraged squared momentum transfer between the medium 
and the fast parton per unit path length q — (q^_)/A. 
This transport coefficient can be calculated from the dif- 
ferential scattering cross section da/d 2 q± in the medium 
or the rate of momentum exchange, dT e i/d 2 q±. In this 
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paper, we will use a common approach to calculate q, us- 
ing rates from Hard Thermal Loop (HTL) effective field 
theory. The basic parameters in the opacity expansion 
are the screening length [ijj and the mean free path A, 
which can also be calculated in HTL. For details, we refer 
the reader to Appendix A 



3. Kinematic approximations in different parton energy 
loss models 

With regard to the technical approximations employed 
in parton energy loss calculations, the four model classes 
show important commonalities, but also display charac- 
teristic differences. 

The three main assumptions that are made in all the 
calculations are: 

• Both the parton and the radiated gluon are on 
eikonal trajectories: the parton energy E is much 
larger than the transverse momentum exchanged 
with the medium q±, £ > q± and the energy of 
the emitted gluon uj is also much larger than the 
exchanged momentum: w > gi. 

• Small-angle (collinear) radiation: the energy of the 
gluon is much larger than its transverse momentum 
fcr, uj ^> kr 

• Discrete scattering centers, or some form of lo- 
calised momentum transfer: the mean free path A 
is much larger than the Debye screening length 1/ \i 
A > l//x. 

Since the eikonal approximation is used, one often uses 
the additional approximation that the gluon energy is 
much smaller than the parton energy w C £, or i = 
uj/E <C 1, i.e. the soft radiation approximation. This 
approximation is not made in the AMY calculation. 

Let us first consider the kinematic constraints E 
uj 3> kx , qi_ , which is often referred to as the soft eikonal 
limit. In the calculations, this limit is used for exam- 
ple to neglect the changes to the parent trajectory due 
to multiple scattering, which simplifies the calculations 
significantly. In phenomenological applications of parton 
energy loss, however, gluon radiation is calculated in the 
entire allowed phase space, up to uj = E and kx = lu, 
where the approximations are not valid. The approxima- 
tions employed in the current calculations lead to finite 
radiation probability at the kinematic bounds u> ~ E 
and kj< uj, thus violating energy-momentum conserva- 
tion. Most current formalisms remedy this by imposing 
explicit cut-offs at the kinematic bounds. The sensitiv- 
ity of the result to large- x and large angle radiation can 
be explored to get an impression of the accuracy of the 
result. The two limits warrant separate discussions. 

Consider the large-x regime, where u> m E. For all 
energy loss models, the gluon energy distribution peaks 
at a "typical" energy (which may be very small). As a 
result, one can always restrict the considerations (and in 



most cases the measurement as well) to a reasonably high 
parton energy for which the calculated rate at x — > 1 is 
small compared to the total rate, thus giving confidence 
that the impact of the cut-off at x = 1 is small. For 
typical medium densities at RHIC, the gluon energy uj 
peaks at 0(1 GeV), so that for parton energies E > 10 
GeV, the probability density at the kinematic boundary 
becomes reasonably small. 

In the present calculational framework, the yield in the 
spectrum for uj > E, i.e. above the kinematical bound- 
ary can be taken as a probability of total absorption of 
the parton, 'death before arrival'. We note that simply 
ignoring the radiation spectrum beyond the kinematic 
limit leads to the unphysical result that in some cases 
the total radiation probability decreases with increasing 
density or path length, when the typical gluon energy is 
close to E. In the AMY formalism, uj > E is kinemati- 
cally allowed because interactions with dynamic medium 
allow the propagating parton to absorb energy from the 
medium when E is not much larger than T, which is then 
re-emitted. 

The situation is different for large-angle radiation. In 
the presently existing formalisms, the typical transverse 
momentum of the radiated gluon k? depends on the typ- 
ical transverse momentum exchanges q± and the number 
of scatterings L/X, but not on the gluon energy uj. As a 
result, there is always some radiation with uj smaller than 
the typical k? and thus with a large probability for radi- 
ation at &t — > oj. The quantitative impact of large angle 
radiation depends on the medium model (large q± con- 
tributions in the medium cross section) and the choice 
of parameters. Already in an early publication it was 
pointed out that the large-angle regime may be impor- 
tant for phenomenological applications, for example in 
|53j : "... the BDMPS-Z formalism is based on the as- 
sumption of small transverse gluon momentum <C uj 
while we find the main contribution to radiative energy 
loss for \kr\ = 0(uj). Both features question the validity 
of the BDMPS-Z formalism ..." A detailed discussion of 
the effect of large angle radiation in the opacity expan- 
sion framework in given in Section III B and in [54] . 



4- Multiple gluon emission 

Multiple gluon emission is calculated by repeating the 
single gluon emission kernel as needed. 

The simplest procedure for multiple gluon emission is 
the Poisson ansatz, where the number of emitted glu- 
ons follows a Poisson distribution, with the mean num- 
ber given by the integral of the gluon emission spectrum. 
The energy distribution of each gluon follows the single 
gluon emission kernel. This procedure is used by GLV 
and ASW. In most cases relevant to the experimental 
conditions at RHIC, the mean number of radiated glu- 
ons (N g ) > 1, so that the mean total energy loss (AE) is 
larger than the mean of the energy of a single gluon (to) . 
In general, this procedure leads to a distribution of lost 
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energy that does not conserve energy as the degrading 
momentum of the parent parton is not dynamically up- 
dated. However, in most cases, the probability to radiate 
a larger energy than the incoming parton energy is lim- 
ited, so that the uncertainty associated with this effect is 
small. For a more detailed study, see [55]. Interference 
effects between medium-induced radiation and vacuum 
radiation are included in the single-gluon emission spec- 
tra, but the parton is assumed to fragment in the vacuum 
after energy loss. 

HT and AMY both use a coupled evolution procedure 
to calculate multiple gluon emission. In the case of HT, 
medium-modified DGLAP evolution is used, which in- 
cludes the virtuality evolution in vacuum. In the AMY 
approach, rate equations are used and no vacuum radia- 
tion is included. However, recent work [56] has addressed 
this issue and showed how it can be included in the same 
framework. 

The evolution equations used in HT and AMY both 
include the coupling between the quark and gluon distri- 
butions in the jets and keep track of the gradual degra- 
dation of the jet energy. The emission probability dis- 
tribution changes as the jet energy degrades, which is in 
principle a significant conceptual improvement over the 
Poisson convolution approach. There is, however, an im- 
portant point that is not explicitly addressed in any of 
the models: as the energy degrades, the remaining path 
length through the medium also decreases. So, in prin- 
ciple, the energy-evolution should be accompanied by an 
evolution in coordinate space and the emission probabil- 
ities should be calculated using local information about 
the parton energy, the medium, and the remaining path 
length. Some models do take into account the evolution 
of the medium and corresponding change in the local en- 
vironment of the propagating parton. However, combin- 
ing both local information and the true finite size effect 
proves to be a difficult problem in the presence of inter- 
ference effects, where the gluon radiation is not a purely 
local phenomenon, but instead couples to the parton over 
an extended area. 



C. Schemes of radiative parton energy loss 
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FIG. 2: The landscape of pQCD based jet quenching for- 
malisms. Arrows indicate common concepts or assumptions 
between adjacent formalisms. 

In the previous sections, we outlined both the single- 



gluon emission calculations as well as the multiple gluon 
calculations that are commonly used in the literature. To 
gain a better understanding of the various formalisms, 
it is useful to compare and contrast closely related for- 
malisms. Fig. [2] illustrates the close technical and con- 
ceptual relations between the different models of parton 
energy loss. 

1. Multiple soft scattering vs Opacity expansion 

The multiple soft scattering approximation and the 
opacity expansion can be shown to be different ap- 
proximations of the same BDMPS-Z path integral. 
In the opacity expansion, the hard scattering tail 
of the scattering potential of the medium is taken 
into account, but the interference between neigh- 
boring scattering centers is only taken into account 
order by order in the expansion. The most com- 
monly used approximation is the single hard scat- 
tering (N = 1) approximation which assumes that 
there is only one dominant hard scattering. The 
multiple-soft scattering approximation resums all 
the interference terms, at the cost of neglecting 
hard scatterings with the medium. In general, the 
multiple-soft scattering approximation is expected 
to be valid for thick media, while the single-hard 
scattering approximation is more accurate for thin 
media. 

2. AMY vs multiple soft scattering approach 
Parton energy loss within the AMY formalism can 
be formulated in terms of the same path integral 
entering the BDMPS-Z formalism [35]. The tech- 
nical commonalities between both approaches are 
further elaborated in 45, 57]. The main concep- 
tual difference between the AMY and the BDMPS- 
Z calculation is the formulation of the medium: an 
equilibrated high-temperature plasma in AMY ver- 
sus static scattering centers in BDMPS-Z. There 
are also important technical differences, however. 
Most notably: in AMY an infinite-length medium 
is used, while BDMPS-Z includes finite-size effects 
in the emission rates. The effect of these differ- 
ences is discussed in detail in [SB]. In addition, the 
ASW calculation include the effect of large angle 
cut-off on the radiation. This effect is not taken 
into account in the standard AMY calculation, but 
is explored in detail in Section |III C 3| below. 

3. Opacity expansions: DGLV vs ASW-SH 

There are two implementations of the (N = 1) 
opacity expansion in the literature: DGLV [33 UBJ 
and ASW-SH [35J[3j5]. Both groups calculate the 
same set of multiple scattering Feynman diagrams 
within the same high-energy and collinear approxi- 
mation. However, in extending this result from the 
kinematic region in which these approximations are 
valid to the entire phase space open for gluon pro- 
duction, they take different approximations. The 
resulting numerical differences are discussed in de- 



tail in Section III B In particular, large angle gluon 
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radiation is a source of significant theoretical uncer- 
tainty. 

4. Higher Twist vs Opacity Expansion 

The Higher Twist formalisms and the (N = I) 
Opacity expansion are related in the sense that 
they are formulated for media where the jet scatters 
a few times per emission. While both formalisms 
have extended the theoretical set-up to include mul- 
tiple scatterings per emission, most comparisons to 
data still make use of the one scattering per emis- 
sion approximation. This is referred to as lead- 
ing order in opacity for GLV scheme and as next 
to leading twist in the higher twist scheme. The 
characterization of the medium is also somewhat 
different as the GLV scheme incorporates the full 
functional form for scattering off a heavy static 
scattering center and the HT scheme takes the 
leading transverse momentum moment of the ex- 
changed transverse momentum distribution. The 
HT scheme thus assumes a Gaussian distribution 
for the exchanged transverse momentum. 

The formalisms differ significantly in the way 
multiplc-gluon emission is treated. The Higher 
Twist calculation uses modified DGLAP virtuality 
evolution, while in the Opacity Expansion, multi- 
ple gluon emission is described using Poisson dis- 
tributed independent emissions (see Section II B 4 1. 
The modified DGLAP evolution used in the HT 
scheme implies that subsequent radiated gluons are 
ordered in transverse momentum fcj^ (kj_ of a given 
radiation is the upper bound for the subsequent ra- 
diation). There is no such restriction in the GLV 
scheme and each radiation may explore the limits 
set by the kinematic bounds. 



D. Summary of the main findings 

A detailed comparison of the different energy loss mod- 
els touches on many aspects of the underlying phsyics as- 
sumptions, but also the technical implementation. Most 
of these aspects are mentioned in the preceding discus- 
sion and are treated in more detail in Section Mil Here 
we list the four most important areas which are funda- 
mental to the problem of radiative energy loss in a QCD 
medium. 

I. Effect of kinematic limits: 

All formalisms make use of collinear approximation 
ui kx- In addition, most of the calculations pre- 



III C). Large angle radiation, however, is expected 



sented in Section III (ASW-SH, ASW-MS, GLV, 



and HT) , except AMY, make use the soft radiation 
approximation E 3> lo. In phcnomcnological cal- 
culations, we have to extend the calculations into 
the large-x, u — > E, and large-angle, kr — > w, do- 
main. For sufficiently large parton energies, the 
uncertainties associated with the large- cc domain 
seem to be reasonably small (see Sections III B and 



to be problematic for all models. This is explored in 
detail in Section [ill B| for the opacity expansions, 
where the resulting uncertainty is found to be about 
a factor 2 in the gluon spectra. The resulting un- 
certainty in Raa decreases with increasing px [54] . 
In Section |III C| it is shown that part of the differ- 
ence between the AMY and ASW-MS formalism 
can be attributed to the absence of a large angle 
kinematic bound in the AMY calculation. 

2. Effect of treatment of multiple gluon emission: 

In the existing formalisms, the calculation of the 
single-gluon emission kernel and the calculation of 
multiple gluon emission are separate steps. GLV, 
ASW-SH and ASW-MS use the Poisson ansatz; 
AMY uses rate equations; and HT uses DGLAP 
evolution. We find that differences between the en- 
ergy loss formalisms originate at the single-gluon 
level. It is possible that the differences between 
the multiple-gluon emission treatment further af- 
fect the total energy loss, but we have not inves- 
tigated this in detail. Semi-inclusive observables, 
like the fragmentation function or the jet shape, 
are expected to be more sensitive to the details of 
multiple gluon emission, such as whether the glu- 
ons are emitted independently (Poisson ansatz) or 
in a kx ordered pattern (DGLAP evolution), than 
the leading particle energy loss or the single hadron 
suppression Raa- 

3. Effect of medium model: 

In most of the current energy loss formalisms, the 
picture of the medium is intrinsically tied to the 
approach used in the calculation. As a result, it is 
not possible to separate the effect of the medium 
model from assumptions made in the rest of the cal- 
culation. For example, the multiple-soft scattering 
calculation ASW-MS is only analytically feasible in 
the limit of gaussian broadening, which corresponds 
to scattering centers with a harmonic oscillator po- 
tential. On the other hand, the AMY approach 
uses a medium based on Hard Thermal Loop field 
theory with dynamical scattering centers. Futher 
work is required to systematically compare differ- 
ent medium models. 

4. Effect of including vacuum radiation: 

All formalisms include vacuum radiation to some 
extent, except AMY. The fact that AMY shows the 
largest suppression in the final comparison (Section 



IV) is mostly due to lack of finite length effects. 
These effects lead to the L 2 dependence of energy 
loss during the early stage of the evolution. At later 
times, or larger L, this turns into a linear depen- 
dence on L. Because AMY implicitly assumes an 
infinite medium the L 2 dependence is lost. This 
behavior is studied in detail in 56J using the light 
cone path integral formalism including finite length 



effects. The obtained rates show the expected L 2 
dependence of the energy loss early in the evolution 
and approach the AMY result at late times. 

Before we end this Section, we emphasize a common 
limitation of all formalisms of parton energy loss dis- 
cussed here: They assume that the hadronization of the 
leading parton occurs in vacuum. The formalisms de- 
scribe the interaction of the leading parton with the 
medium strictly at the partonic level, but do not con- 
sider any influence of the medium on the hadronization 
process or final-state interactions of the produced leading 
hadron with the medium. This appears to be justifiable 
for leading hadrons with momentum ph S> rrih, where 
to/j is the hadron mass, because the formation time for 
such a hadron is Lorentz dilated by a factor Ph/rrih ^> 1, 
and thus the formation of the hadron may be assumed 
to occur far outside the medium. Indeed, as discussed 
in the introduction, experimental evidence supports this 
assumption for light hadrons. We note however that this 
assumption is not generally satisfied under present con- 
ditions of heavy ion collisions at RHIC for baryons and 
mesons containing heavy quarks. 



E. Directions for future work 

It is clear from the preceding discussion that the cur- 
rent energy loss calculations have a number of specific 
weaknesses which may have large quantitative impact on 
the results. Based on our review of the existing models, 
we can formulate a number of questions for future work. 

Firstly, we recommend to develop calculations with 
better control of the large angle and large energy radi- 
ation. In a realistic calculation, the radiation cross sec- 
tion should naturally go to zero at the kinematic limits 
kx = lu and u) = E. It is most important to control 
the large-angle behaviour; the large-energy regime can to 
some extent be avoided by concentrating on high-energy 
partons/jets. The most natural way to control the large- 
angle behaviour of the radiation would be to use the full 
Next-to-Leading Order (NLO) matrix elements for 2 — > 3 
scattering with radiation. 

The treatments for multiple gluon radiation should 
also be further investigated. So far, the treatments are 
based on incoherent superposition and/or fc^-ordering of 
the shower. It is worthwhile to quantitatively investigate 
deviations from these assumptions. Semi-inclusive ob- 
servables such as fragmentation functions and jet shapes 
may be more suitable to investigate specific multi-gluon 
emission scenarios than inclusive observables, such as the 
nuclear modification factor Raa- 



III. DETAILED COMPARISONS OF MODELS 

In the following sections, we will present more detailed 
side-by-side comparisons of the existing energy loss for- 



malisms. The discussion follows the structure layed out 
in Section HTCl 



A. The ASW formalism 

As discussed in Section [II B 2| the ASW formalism cal- 
culates parton energy loss based on a path-integral for- 
malism 38 - 40 which can be evaluated either in the mul- 
tiple soft scattering limit or through an opacity expan- 
sion. The GLV N — 1 opacity result reproduces the 
ASW expression [35J on the level of the Feynman dia- 
grams and the analytic expression for the ui- and kx- 
differential gluon energy distribution. 

In this section, we concentrate on the multiple soft 
scattering approximation ASW-MS, which has been 
mostly used in comparisons to experimental data. A 
comparison of the single-hard scattering result ASW-SH 
to the DGLV opacity expansion g5] used by WHDG [27] 
is given in Section [HI B| 



1. Gluon radiation spectrum 

The result of the ASW-MS calculation is publicly 
available in terms of the 'quenching weights' which rep- 
resent the energy loss probability distribution [58 . How- 
ever, before discussing the probability distribution, we 
will present the single gluon radiation spectra. The mul- 
tiple soft scattering limit suppresses the production of 
infrared gluons by a destructive interference effect. As a 
consequence, all spectra are peaked at finite gluon ener- 
gies. In general, the radiated gluons become harder as 
one increases the average energy loss (i.e. as one increases 
q). If the projectile energy is sufficiently large and the 
in-medium pathlength is sufficiently small, then the radi- 
ated gluons carry small fractions of the projectile energy. 
This is the case for a projectile quark energy E = 100 
GeV shown in the upper panel of Fig. [3] 

However, if the projectile energy is too small, see lower 
panel of Fig. [3] one faces a particular problem: One 
calculates the radiated gluon spectrum as if the parton 
would propagate through a medium of path-length L, al- 
though there is a finite probability that the parton does 
not have sufficient initial energy to make it through L, 
and gets stuck before. In the present calculational frame- 
work, finding yield in the spectrum beyond the kinematic 
boundary x = 1 (u> = E), signals that one has assumed 
that the particle propagates through a length L for which 
the probability of 'death before arrival' is finite. 



2. Multiple gluon emission 

In the small x limit, the gluon emission probability 
is independent of the incident parton energy. In this 
case one may convolute with the Poisson distribution to 
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FIG. 3: Energy spectrum of radiated gluons, for a light quark with energy E — 100 GeV and path length L — 5 fm (top panel) 
and for a light quark with E = 10 GeV and L — 2 fm (lower panel). The legends on the plots indicate the average energy loss 
and the corresponding value of the transport coefficient q. 
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FIG. 4: Energy loss probability density P(AE) for a light quark of energy E — 100 GeV and path length L = 5 fm (top panel) 
and E — 10 GeV and L — 2 fm. The different lines are for different average energy loss, as indicated in the legend, together 
with the corresponding values of the transport coefficient q. 



compute the probability P(AE) for an incident parton to 
lose energy AE during its passage through the medium 



71=0 ' 1=1 



duj 



dl 

dcu 



S(AE-J2^),(1) 



i=l 



where (N g ) — J °° dui(dl / duj) is the mean number of ra- 
diated gluons. This Poisson convolution is used in the 
ASW calculations for the multiple soft gluon approxima- 
tion as well as the opacity expansion, and also in the 
GLV-based opacity expansions. 

The resulting energy loss probability distribution 



10 



P(AE/E) are shown in Fig. [I] using the ASW multiple 
soft scattering limit for an incoming quark with energy 
E and in-medium path length L. All results have been 
computed with a s = 0.3. The probability distribution 
has three distinct pieces: 

1. Untouched survival: There is finite probability that 
a parton does not interact with a medium of length 
L and that it loses no energy. This probability is 
P(0) = exp(— (N g )) and is represented by a color- 
coded dot at AE/E = -0.05 in Fig. [I] 

2. Survival with finite energy loss: This is the con- 
tinuous probability that a parton makes it through 
a medium of length L but loses a finite fraction 
AE/E during its passage. This is denoted by the 
color-coded curve at finite < AE/E < 1. 

3. Death before arrival: In general, if one shoots a par- 
ticle into a wall of thickness L, it can get stopped 
on its journey before reaching the length L. This 
probability is denoted by the color-coded dot at 
AE/E = 1.05. 

In general, as the average energy loss grows due to an 
increase in q, one observes that: 

(a) the probability of untouched survival decreases; 

(b) the probability of survival with finite energy loss 
shifts to larger values of AE/E; 

(c) the probability of death before arrival increases. 

This is seen clearly in Fig. [4j which shows the energy 
loss probability distribution P(AE/E) for two different 
choices of the parton energy E and path length L. The 
extreme case is E = 10 GeV and L = 2 fm (lower panel 
of Fig. [1]), where requiring an energy loss of AE = 4 
GeV amounts to a greater than 40 percent probability of 
untouched survival and 30 percent probability of death 
before arrival. This comes close to an all-or-nothing 
scenario, where a particle either goes through without 
medium modification or gets stuck, but its probability of 
emergence as an object with reduced energy is relatively 
small. 



3. Uncertainties associated with kinematic limits 

The limitations of the high-energy eikonal approxima- 
tions used in the derivation of the path-integral formal- 
ism from which both the multiple soft scattering limit 
and the opacity expansion stem, were discussed in the 
original papers [58H4T)] . More recently, ASW have used a 
number different approaches which are trying to address 
the shortcomings of the analytical calculations, such as 
using DGLAP evolution equations [35] and Monte Carlo 
algorithms for final state radiation [BDll63| . 




o.o 0.2 0.4 0.6 0.0 1.0 



x 

FIG. 5: Single inclusive gluon radiation distribution, 
dN g /dx, from the WHDG implementation of the first order 
in opacity DGLV formula, Eq. S, in red, and the ASW-SH 
implementation of Eq. {3j, in black, for a 10 GeV up quark 
traversing a nominal, 2 fm long static brick of QGP held at 
a constant T = 485 MeV. The point at x = 1 indicates the 
integrated weight of dN g /dx in the ASW-SH implementation 
for x > 1. 



B. Opacity expansion: WHDG and ASW-SH 

This subsection contains a detailed comparison of the 
radiative part of the WHDG calculation [47] . equivalent 
to the first term in the DGLV opacity expansion |46j . 
and the ASW opacity expansion formalism truncated to 
first order in opacity, here called ASW-SH [3H ED] . We 
discuss the importance of the definition of the splitting 
variable x and of the kinematic cut-offs made to enforce 
the assumptions used in deriving the energy loss formu- 
lae. Quantitative results comparing the radiative com- 
ponents of WHDG and ASW-SH will also be shown. 



1. Single emission kernel 

One of the major driving forces of the present investi- 
gation was the realization that some of the approxima- 
tions made to arrive at analytical results for the medium- 
induced radiation have a larger impact on the final re- 
sult than expected. An excellent case in point comes 
when one attempts a naive comparison between the sin- 
gle inclusive gluon distribution implemented in a massless 
quark and gluon version of the radiative piece of WHDG 
and the one found from the ASW-SH code pM [55] for 
the opacity expansion. Both purport to compute the sin- 
gle inclusive distribution of gluon radiation, dN g /dx, to 
first order in opacity [55J HU [52] for a medium of Debeye- 
screened colored static scattering centers [TTj . See Fig. [5] 
which shows dNf GlN /dx and dA^ sw_SH /dx for a nom- 
inal 10 GeV quark jet in a static plasma of length 2 fm 
and temperature T = 485 MeV [HI]. 
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The DGLV formula for the first order in opacity energy loss is 



dN T>GLV 

dx 



2C R a s L 

7T 3 A 



d 2 qd 2 k- 



kq(k-q) 2 - P 2 (\- (k-q) 



2 ^ [(k-q) 2 + /3 2 ] 2 (k 2 



P 2 



dz 



cos 



(k-q) 



2xE 



p{z). 



Here q is the transverse momentum exchanged with the medium and k the transverse momentum of the radiated 
gluon, fi is the Debeye screening mass, m g and M q are the effective thermal mass of the radiated gluon and the mass 



of the parent parton, respectively, j3 = m 



~M 2 ; A is the mean free path of the radiated gluon and p(z) is the 



probability distribution for the distance to the first scattering center 
which docs not include the effect of a thermal gluon mass m g , is [40] 



j. The equivalent expression for ASW-SH, 



^jASW-SH 



dui 



2Cna s L 

7T 3 A 



d 2 qd 2 k 



dz 



P 



k • q(k 



x 2 M 2 q • (k - q) 



(q 2 



[(k - q) 2 + x 2 M 2 ] (k 2 



cos 



(k- 



x 2 M 2 



2lu 



x 2 Ml) 



(3) 



If we set u = xE in (J3j) and m g — in j3 in 
then we see that Eqs. ^ and ^ are, in fact, identical. 
Why, then, do the curves obtained by evaluating these 
expressions, shown in Fig. [5] differ so drastically? The 
answer lies in the details of the implementation of these 
two expressions. It turns out that there are a number of 
choices that must be made in going from the differential 
expressions for the radiation spectrum to the integrated 



gluon distribution, 
Table H 



dNg/dx. These are summarized in 





WHDG 


ASW-SH 


^max 


2x{\ - x)E 


xE 


X 


x+ 




*7max 


y/3pE 


oo 


Pi*) 


2e- 2z/L 9(z)/L 


6(L-z)0(z)/L 


L/X 


Lpo = Lp x 9-Kaj/(2p, 2 ) 


1 


a s 


0.3 


1/3 


m g ; M q 


H/y/2; n/2 


0; 



TABLE I: Table of differences between the implementations 
of Eqs. (2J and ^ in WHDG and in ASW-SH. 

The differences outlined in Table U have been individ- 
ually explored to see which choices have a large effect 
on the calculated gluon spectra. The effect of non-zero 
quark mass M q is small, because the mass is small com- 
pared to the total energy. The choices for as differ by 
10% which directly translates into a change in the single 
gluon spectra by the same amount. The shape of the 
distribution also depends on as through the value of p, 
but this is a small effect. It also turns out that the choice 
of Qmax has only a small effect, because of the small radi- 
ation probability for q > y/3fiT. Another small effect is 
the choice of p(z); the difference between the exponential 
profile and a uniform density only leads to characteristic 
shapes in the spectrum for small L < 2 fm, but does not 
effect the total radiation probability much. The effect of 



T = 485 MeV, L ; 
E = 10GeV 



2 fm 



2x(1-x)E, m , m 



k max =xE, m =m^=0 
k m ax =x E, L/l=1, 




FIG. 6: Comparison of the single gluon emission spectrum 
for several variants of the N = 1 opacity expansion, using 
a uniform medium with T = 485 MeV and L = 2 fm. All 
curves use as = 0.3 and g max = 10 GeV 3> \/3fiE. The 
solid curve labeled ASW corresponds to what is used in [391 
140] . while the grey (red) curve corresponds to the radiative 
WHDG formalism used in [461 Wf\ , except for the difference 
in (/max- See text for details. 



the density profile p(z) on the gluon spectrum can be seen 
in Fig. |6j comparing the lower solid line, labeled ASW, 
which uses p(z) = 9(L — z)/L and the dot-dot-dashed 
line using p(z) — 2exp(—2z/L)/L. 

The remaining three items in Table [I] have signifi- 
cant effects on the gluon radiation spectrum. The ef- 
fect of each of these on the gluon spectrum is show in 
Fig. [6j L/X enters as a multiplication factor in the sin- 
gle gluon spectrum (compare dot-dot-dashed and dashed 
line). This is a significant factor at T = 485 MeV, 
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A = 0.62 fm. The choice of L/X = 1 in the ASW papers 
|40j was made for convenience in tabulating the quench- 
ing weights P(AE), but this calculation has never been 
used to determine the medium density. In the rest of this 
paper, we will therefore include L/X different from unity 
in the calculation. Introduction of a finite gluon mass m g 
leads to a sizable reduction of the radiation at small x 
where uj w m g (compare dashed and dotted lines in Fig. 
[6]). On the other hand, due to nontrivial interference ef- 
fects, its effect on observables such as Raa is actually 
not very large [M] , 

Finally, the choice of the cut-off k max on the perpen- 
dicular component of the gluon momentum has a large 
effect, mostly at small x, where the difference between 
the cut-offs is a factor 2 and the probability of radiation 
at finite k± is large (compare dotted line and solid grey 
(red) line in Fig. [6]). Because this cut-off is also related 
to large angle radiation, it has a fundamental impact on 
our ability to calculate the radiation spectrum, which we 
will discuss in more detail in the following section. 



cut-off angle # max can be defined, leading to 

, _ J x + £; + tan(6» max /2), x = x+, 
y x E Esm{0 ma ,x), x = x E . 

Note that the same physical cut-off at an angle 9 max 
gives different expression for fc max as a function of xe 
and x+. The exact relations between x+ and x E are 

Note that to lowest order in the expansion parameter 
is Ut/xE, which controls the extent of collinearity, the 
two definitions of x are identical. The difference between 
using x e and x + in the calculations is only significant at 
larger angles. 



2. Definition of the variable x 

The upper bound of the transverse gluon momentum 
k± is different in the WHDG calculation than in the 
ASW-SH calculation. At small x, both formalisms im- 
plement a cut to ensure forward emission of the gluon. 
The difference in the values of the cut-off arises be- 
cause in the derivations descending from the work by 
GLV[HJ EH IH x is taken to be the positive light- 
cone momentum fraction x = x + = k + / P + , while the 
ASW-SH derivation OED] uses x = x E = k°/P°, the 
fraction of energy carried away by the radiated gluon. 

If we denote the usual four-momenta with parentheses 
and light-cone momenta with brackets then the radiated 
gluon four-momentum is 



k = (x E E,k, V(x E £) 2 -k 2 
k 2 



x+E+, 



x+E^ 



(4) 



where k is the momentum of the gluon transverse to the 
direction of the parent parton and the gluon is assumed 
to be on mass shell. 

In light-cone coordinates forward emission implies that 
k + > kr . This condition implies 



kx < k„ 



.1" 



(5) 



For a massless parent parton, E + = 2E, so that fc max = 
2x+E. In Minkowski coordinates forward emission im- 
plies that k z > 0. This condition implies 



fcr < k n 



x E E. 



(6) 



3. Transverse momentum cut-off at large momentum 
fraction 

Eqs. ([2| and Q also have support for all values of x. 
Nonzero weight in dN g /dx for x > 1 of course violates 
energy-momentum conservation. Requiring the contin- 
ued forward propagation of the parent parton leads to 
an additional kj- cut-off that enforces energy-momentum 
conservation and also restricts large-angle radiation at 
large x. 

The momentum of a massless parent parton after en- 
ergy loss is 



(l-x E )S,q-k, V / ((l 
(q-k) 2 



x E )E) 



(q-k) 



(i 



(1 - x+)E^ 



(10) 



where q is the transverse momentum transfer to the par- 
ent parton from the scattering center. For completeness 
the original parent parton momentum is P — (E, 0, E) = 
[£7+0,0]. 

In light-cone coordinates, forward propagation of the 
final state parton implies p + > p~; in Minkowski coor- 
dinates it implies k z > 0. For the light-cone coordinate 
case forward propagation leads to 



(1 



)E+ >|q-k| 



(11) 



where |q| ~ 3T < q max = V6ET is small compared to 
most values of |k|; for the Minkowski coordinate case 
forward emission leads to 



(1 - x E )E > |q - k| 



(12) 



Both values of fc max restrict emission to angles be less For each definition of x there are two cut-offs (e. g. 



than 90°, which is still a rather wide angle. In general, a Eqs. ^ and (11 1 for 



or Eqs. ^ and Q for 
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x + 

FIG. 7: Single gluon spectra using Eq. S with fc max = 
x+E + (black) and fc max = x+(l — x+)E + (red) cut-offs for a 10 
GeV up quark traversing a 2 fm static QGP of T = 485 MeV 
and m q = m g = 0. The black dot at x+ = 1 represents the 



integrated weight of dN g /dx+ for x+ > 1 when k n 



■■ x+ 



E+. 



x = ie); this needs to be taken into account when eval- 
uating the kr integral over dNg/dxdkx- One possibility 
would be to take the minimum of the two definitions; 
for instance, using light-cone coordinates and taking 
#max = tt/2; this would mean fc max = min(a; + , 1 — x+)E + . 
The present implementations of DGLV and WHDG use 
a smoother function, fc max = x + (l — x + )E + . Note that 
the existing ASW-SH implementation [55] does not 
include a large-a; cut-off. 

Figure [7] shows the single gluon spectra dN/ dx g with 
the two different cut-offs Eq. [5] and Eq. [IT] It can be seen 
in the Figure that the effect of introducing the kinematic 
cut-off x < 1 is small in regimes where the dN g /dx peak 
is large and at small x, for instance for L = 2 fm, with 
the effective gluon mass set to zero, m g = 0. In this 
sense the dN g /dxdkrdqx integrand respects the small- a; 
approximation rather well. In Fig. [6] the effect of lifting 
the x < 1 kinematic cut-off is more important because 
the number of radiated gluons (N g ) is smaller when m g ^ 
0; in this case the large- x cut-off causes a 10-20% change 
in (N g ). In general, the effect of the large-a; cut-off is 
expected to increase with increasing path length and to 
decrease with increasing parton energy. 



4- Effect of large angle radiation 

In order to illustrate the sensitivity of the single gluon 
emission probability to fc max , Fig. [8] plots dN g j dxdkr for 
x = 0.025, along with an illustration of three possible 
cut-offs for kr'- 

(a) x — x+ with 6 max = 7r/2; 

(b) x = x E with max = 7r/2; 

(c) x = x E with max = 7r/4. 

Recall that to lowest order in collinearity, the first two 
cut-offs are identical; the third is a natural 0(1) variation 




FIG. 8: Transverse momentum spectrum dN g /dxdk,T of 
emitted gluons with x — 0.025, calculated using Eq. Q for a 
light quark with all masses set to 0, E = 10 GeV, L = 5 fm, 
and representative values of fi w 0.46 GeV and A ~ 1.25 fm 
for a medium density of dN g /dy = 1000 similar to RHIC con- 
ditions [47] . Vertical lines depict the three values of kr dis- 
cussed in the text as possible cut-offs to enforce collinearity 
in Eq. Q. 



in the cut-off that one can use to estimate the systematic 
uncertainties deriving from the collinear approximation. 
Clearly the assumption of collinearity is badly violated: 
for values of x ~ fJ-/E, dN g /dxdkT reaches its maxi- 
mum value at kr ~ xE. For these values of x the emis- 
sion spectrum is highly sensitive to the choice of fc max : 
dN g /dx - *4 ax . 

Since the collinear approximation is so badly broken, 
it is not a good approximation to take 5G_|_ k ig. A 
meaningful comparison of results, then, can come only 
when the emission spectra of Eqs. ^ and ^ are plotted 
with respect to the same variables. Since one is interested 
in a differential quantity, a Jacobian is required. We 
choose to transform x+ to ie because, ultimately, one is 
interested in energy loss, as opposed to the loss of positive 
light-cone momentum. The transformed spectrum is then 
given by 



dN L v 

(ze) 



dxE 

with 



dki 



dx + dN g 
dxE dx + dkr 



(x+(x E )), (13) 



dx+ 
dxE 



1 

2 

= xeE sin( 



\x E EJ 



(14) 
(15) 



Note the change in the upper limit ( 15 ) of integration in 



Eq. jl3|). The resulting comparison of dNg/dx^ is shown 
in Fig. |9| Note the very large difference in the results for 
the two collinearly equivalent definitions of x and that for 
the result with a reduced 6* max . Of course this enormous 
difference implies a large (factor 2 — 3) uncertainty in 
the extraction of the medium parameters from leading 
hadron suppression data |54j . 
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FIG. 9: Comparison of Eqs. |2| and Q in the massless 
limit and for which the x+ dependence of Eq. |2| has been 
transformed into ie; see Eq. ( |13[ ). Also shown is the result 
when using the Xe interpretation and reducing # max to n/4, 
a reasonable 0(1) variation in the kr cut-off. 



C. AMY, BDMPS-Z and ASW-MS 

1. AMY transport equations 

The medium-induced radiative energy loss suffered by 
high energy partons passing through nuclear matter was 
first computed in the BDMPS-Z approach [H I54U57], 
in which the gluon emission probability is expressed in 
terms of the Green's function of a 2-D Schrodinger equa- 
tion with an imaginary potential proportional to the 
interaction cross section with color centers of a quark- 
antiquark-gluon system. 

In the AMY approach [51, 52, 65, 66 , the gluon emis- 
sion rates are calculated fully at leading order in a s by 
resumming an infinite number of ladder diagrams in the 
context of hard thermal loop resummed QCD. Both ap- 
proaches are valid in the multiple soft scattering limit, 
but differ in several essential ways: in AMY the medium 
consists of fully dynamic thermal quarks and gluons, 
while in BDMPS-Z the medium is treated as a collec- 
tion of static scattering centers. In BDMPS-Z the gluon 
emission probability is calculated in configuration space 
while in AMY the radiation rate is calculated in mo- 
mentum space. Salgado and Wiedemann [39 further ex- 
tended the BDMPS-Z formalism to include the correct 
thin plasma limit, which has important quantitative ef- 
fects. In addition, different evolution schemes are used 
for multiple-gluon emission: the AMY formalism uses 
rate equations to obtain the final parton distributions 
while applications of the BDMPS-Z calculation convo- 
lute the radiation rate with a Poisson distribution to ob- 
tain the quenching weights. 

The main assumption in these two formalisms is that 
the temperature of the medium is high enough such that 
the asymptotic freedom of QCD makes it possible to treat 
the interactions between a fast parton and the medium 
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FIG. 10: A typical diagram calculated in the AMY and 
BDMPS-Z approaches. 

using perturbation theory. In this case, soft exchanges 
between the medium and the propagating parton domi- 
nate the stimulated radiation of a hard collinear gluon. 
At the same time, the effect of multiple collisions is re- 
duced due to the coherence between multiple soft scat- 
terings within the formation time of the emitted gluon 
(LPM effect). This effect makes it necessary to resum all 
diagrams as depicted in Fig. [I0]to calculate the leading 
order gluon emission probability/rate. 

In the AMY approach, one considers a hard parton 
traversing an extended medium in thermal equilibrium 
with asymptotically high temperature T — > oo. Due to 
the small coupling g — > 0, a hierarchy of parametrically 
separated scales T > gT > g 2 T makes it possible to con- 
struct an effective field theory of soft modes (modes with 
momentum |k| ~ gT) by summing contributions from 
hard thermal loops into effective propagators and vertices 
[BTl . The hard parton traversing a thermal QGP un- 
dergoes a series of soft elastic scatterings with transverse 
momentum of order ~ gT off the thermal particles of the 
medium. The differential cross section (interaction rate) 
at leading order in a s is 

dr c i i g 2 Tml 

d 2 q ± (27r)2qi(qi+ma)- 1 ° J 

Note that the rate has been divided by the quadratic 
color Casimir Cr of the parent parton, indicated by plac- 
ing a bar over the rate r e i (similarly for other quantities) . 

These soft multiple scatterings induce collinear split- 
ting of partons. The time scale over which the parton 
and emitted gluon overlap is of order y/ui/q, which is 
of greater or equal order of magnitude than the mean 
free time of soft scatterings for lj > T, with u> the en- 
ergy of the radiated gluon. To obtain the leading-order 
gluon emission rates, one must consistently take into ac- 
count the multiple scattering processes. Within the ther- 
mal field theory, one essentially calculates the imaginary 
parts of an infinite number of gluon self-energy ladder di- 
agrams. The resummation of these ladder diagrams can 
be organized into a Schwinger-Dyson type equation for 
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FIG. 11: Diagrammatic representation of the Schwinger-Dyson equation for the emission vertices. 



the dressed radiation vertex, as depicted in Fig. [TT] The 
corresponding integral equation is 



2h = iSE(h,p,k)F(h) + J d 2 q x 



d 2 q± 

(C s -C A /2)[F(h)-F(h-fcq ± ) 
-(C A /2)[F(h)-F(h+pq ± )] 
-(C A /2)[F(h)-F(h-(p-A)qi)] . (17) 



In the above equation, 5E is the energy difference be- 
tween the initial and final states. 



SE(h,p,k) = Ep- k + E k — E, 
h 2 m] 



(18) 



l p-k 



2pk{p-k) 2k 2{p-k) 2p ' 



with p the incoming parton momentum and k the mo- 
mentum of the radiated hard gluon. 1/5E is of order 
of the formation time for the bremsstrahlung process in 
the medium. The masses appearing in the above equa- 
tions are the medium induced thermal masses. The 2- 
D vector h is a measure of non-collinearity of the final 
states, defined to be h = (p x k) x e||, with e|| the unit 
vector along the chosen longitudinal direction; here h is 
parametrically of order gTE{w /T) 1 / 4 and therefore small 
compared to p • k, with uj ~ k. For the case of g — > qq, 
the (C s — Ca/2) term is the one with F(h— pq±) rather 
than F(h— fcqj.) in the above equation. 

The gluon emission rate dT{p,k) / dk is obtained by 
closing off the dressed vertex with the bare vertex in- 
cluding the appropriate statistical factors, 



dT(p, k) = C s g 2 1 1 

dk Wnp 7 l±e- fe / T l±e-b- fe )/ T 



l + (l-x) 2 
X 3 (1 — x) 2 



N { 



x 2 + (l-xf 
X 2 (1 — x) 2 
l+x 4 +(l-x) 4 
x 3 (l-x)3 



q q<j 
g -> qq 
g -> 99 



> x 



d 2 h 



2h-ReF(h,p,fc). (19) 



Here C s is the quadratic Casimir relevant for the process (Cf = 4/3 for quarks or Ca = 3 for gluons), and x = k/p 
is the momentum fraction of the gluon (or the quark, for the case g — > qq). The subsequent multiple emissions 
are treated by evolving the momentum distribution P(p) = dN/dp of the traversing hard partons in a set of rate 
equations, 



dt 

dP g {p) 

dt 



= / P q<j {p+k) 



P q q{y+k) 



dT* g (p+k,k) 

dk 

dT$ g (p+k,p) 
dk 



dr q %( P ,k) 



dk 

dT/ g (p+fc,fc) 
dk 



2Pg( P +k) 



d^qq(P+k, k) 



dk 



P 3 (P) { 



(dT 9 qq (p,k) drf g (p,k) 



dk 



dk 



(20) 



Q{k-p/2) . (21) 



In the AMY approach, the medium consists of fully dynamic thermal quarks and gluons; the effect of emitting and 
absorbing thermal energy are fully included in the evolution equations. In addition, the elastic energy loss due to the 
recoil of the dynamical scattering centers can also be consistently included in the formalism p9l I70j . However, the 
transition rates are calculated in momentum space assuming the thermodynamic limit, i.e., the high energy parton 
experiences a uniform medium on the time scale of the formation time of the emitted radiation. 



2. Comparison of BDMPS-Z and AMY 

In the BDMPS-Z formalism, one calculates the amplitude for a quark-antiquark-gluon system evolving in the 
medium without inducing inelastic reactions (following Zakharov). The general formula for the probability of the 
gluon emission by a hard parton traversing the medium may be expressed as, 

dl nrP _^ (t) f°° f°° / \ 

r n li2 Re / #i / *2(V bl -V b3 [G(b 2) t 2 |b 1 ,t 1 )-G vac (b 2 ,i 2 |b 1) t 1 )]) . (22) 

UK [X(i~X)p\ Jq J t \ /bi=b 2 =0 



16 



Here I is the probability of gluon bremsstrahlung from 
the high energy particle, P s ^ g {x) is the vacuum split- 
ting function for relevant process and G(b 2 , i 2 |bi, t-y) is 
the Green's function of a 2-D Schrodinger equation with 
Hamiltonian 

H(p h ,b,t) = SE(p b ) - iT 3 {b,t) . (23) 

The initial condition for the Green's function is 



G(b 2 ,i|b 1 ,i) = 5 2 (b 2 -b 1 ). 



(24) 



The kinetic term in the Hamiltonian describes the energy 
difference between initial and final states, as given by 
Eq.(18l, with pb = h/p, and T3 in the potential term is 
the 3-body interaction rate, 



r 3 (b,i) = ^C A T 2 (b,t) + (C s -^C A )T 2 (xb,t) 



1 



c A r 2 ((i - x)b,t) 



(25) 



where r 2 is related to the Fourier transform of the elastic 
collision rate dT c \/d 2 q± by 



F 2 (b,t) - / d 2 q ± 



1 



(26) 



Note that in the original BDMPS-Z formula, the rate f e i 
for soft scatterings was written as the number density 
p of the static scattering centers in the medium times 
elastic cross section a c \ for such scatterings. Here we 
follow the notation of Arnold [57] and write BDMPS- 
Z formulas more generally in terms of the rate r e ] for 
elastic scatterings off the medium constituents. Since 
the medium is treated as a collection of static scattering 
centers, one replaces the denominator of Eq. (16) with 

(qi + < 



for the elastic scattering rate. Further as- 
suming that the gluon emission rate is dominated by the 
region b < l/mo, one may approximate the elastic rate 
by 



f 2 (b,i) = ^ 2 . 



(27) 



where Cjiq is transverse momentum squared transferred 
to incident parton per unit time (Cr is the appropri- 
ate color factor for the parton). Neglecting the effective 
masses of the particles, the Hamiltonian takes the Har- 
monic oscillator form, 



where 



.(l-x)C A + x 2 C s 



2M 



(28) 



(29) 



and M = x(l — x)p. Making use of the oscillator Green's 
function and performing the integration over time t± and 



t 2 for a brick of medium with length L, one may obtain 
the BDMPS-Z formula for multiple-soft scattering, 



dl axPg^gix) 
dk 7r 



In | cos(wq-^)| 



(30) 



The correspondence between Eq. 22 which is the 
starting point to derive the BDMPS-Z formula above, 
and Eq. [19] the rates calculated by AMY, has been ex- 
plored by Aurenche and Zakharov for photon radiation 
[TTj and by Arnold for gluon radiation [57]. Here we 
reproduce the main steps. First, note that the Green's 
function only depends on the time difference At = t 2 — 1\ 
due to time invariance. Performing the integral over t\ 
which gives a factor of total time, the resulting non- 
vacuum part of Eq. ( 22 ) becomes 



(31) 



dT _ aP s ^ g {x) 
dk x 2 (l — x) 2 p 3 

x Re / dAt [V bl ■ V ba G(b 3 ,Af|bi,0)] b b 2=0 
Jo 

One may define the time-integrated amplitude 

f(b) = 2i/ dt[V bl G(b,t|bi,0)] bl =o, 
Jo 

which satisfies the following equation, 
-2V b <5 2 (b) = Hf(b). 



(32) 



(33) 



In terms of the amplitude f(b), the gluon emission rate 
becomes 



dT 

dk 



aP s ^ g {x) 
x 2 (l — x) 2 p 3 



Re[(2i)- 1 V b • f(b)] b= 



o ■ 



(34) 



Fourier transforming Eq. ( 33 1 from impact parameter 



space to momentum space with the use of Eq. ( 28 ) , one 
obtains the linear integral equation [Eq[l7], with F(h) = 
pf(p b ). The gluon emission rate becomes 



dT _ aP s ^ g (x) 
dk ~ 4x 2 (l - x) 2 p 7 



d 2 h 



Re[2h • F(h)] . (35) 



Including the appropriate final state statistical factors, 
one reproduces the AMY formula ( 19 ) for the gluon emis- 
sion rate. 



3. Large angle radiation in AMY 

In order to explicitly exclude large angle radiation, 
which is automatically suppressed for g <C 1 but can 
contribute for larger g, we suppress the large angle con- 
tributions in the integral over h in Eq. [l9]by introducing 
an envelope function W(k±, fcj_ max ), where k± = h/p. 

A reasonable choice for W is a Gaussian that cuts off 
kr above a value k± max . We choose a Gaussian because 
the integral equation Eq. [17] is solved by Fourier trans- 
forming into the so called impact parameter space (see 
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FIG. 12: (Color online) Envelope function W for different 
maximal transverse momenta. 



[72]). where it becomes a differential equation. Hence, 
we want an envelope function whose Fourier transform is 
analytic and well behaved. Explicitly wc dchnc 



W(k±,k ±max ) = cxp(-kj_ Ai max ) • 



where we can vary k± max , depending on how much we 
want to restrict the emission angle. First, we would like 
to exclude the k± that are kinematically disallowed, so 
we choose the following k± max : 




> k 
< k 



for \p\ > |fc| and \p — k 
for \p\ > \k\ and \p — k t 
for \p\ < \k\ and \p\ > \p — k\ 
for \p\ < \k\ and \p\ < \p — k\ 

C 



(37) 



This suppresses most of the kinematically disallowed 
kj_ and makes sure that k± is always smaller than the 
smallest momentum. Note again, that in the limit of g <C 
1 those kinematically disallowed transverse momenta are 
already suppressed, so that the formalism is perfectly 
consistent with kinematics in the limit where it is derived. 

To also suppress k± which correspond to emission an- 
gles larger than a certain 6* max , we can replace fc_Lmax by 
fcimax sinf^max)- Fig. 12 shows the envelope function W 
for the original fcj_ m ax and for a maximal emission angle 
of #max = 7r/4. Using a Gaussian, not all contributions 
from k± > /cj_ m ax are removed but they are sufficiently 
suppressed to give a good estimate of the effect of remov- 
ing these contributions. 



4- Modified rates 

We demonstrate the effect of suppressing large trans- 
verse momentum emissions by showing the ratio of the 
rate including the suppressing envelope and the original 
AMY rate in Fig. [13] as a function of k, the absolute 
value of the three-momentum of the emitted gluon. 

For a typical a s of 0.3, the suppression at small k is 
significant but decreases for larger k. As expected, the 
effect is also reduced as a s decreases, as can be seen in 
the result for a, = 0.01. 



5. Comparison AMY, BDMPS-Z, ASW-MS in a brick 

Figure [M] compares the gluon emission spectra for 
AMY, BDMPS-Z and ASW-MS for two path lengths 
(X = 2 and 5 fm) and T = 300 MeV. In the AMY formal- 
ism the transport coefficient q, which is used to set the 
medium density in the BDMPS-Z and ASW-MS calcula- 
tions, is a derived quantity not a parameter to be set. For 
this comparison q was calculated usi ng the q T limit 
of the HTL scattering rate (Eq. (16), with <j, 2 nax = ET, 
so that q = 6a s T[i 2 ET / p?))) . 



The AMY curve (green line in Fig 14 1 is obtained 
by multiplying the rate dT/dk Eq. 19 by the medium 
length L to obtain a radiation spectrum dN/dk. The 
curve labeled BDMPS-Z (black dashed) was calculated 
~ while the 'BDMPS-Z (low-x)' curve (solid 



(36) us ing Eq. [30 



black) was obtained using in addition the small- a; approx- 
imation xP q ^ q (x) rs 2 and loq = —iCAq/zp- Comparing 
first the AMY curves with the BDMPS-Z curve in Fig. 
[14} one sees that the results are similar for L = 5 fm, 
while for short path lengths (X = 2 fm), AMY gener- 
ates more radiation. This is due to the fact that AMY 
uses the limit L — > oo which ignores the increased ef- 
fect of finite formation times at small length [73 . The 
ASW-MS calculation (blue curve in Fig. 14 ) is based on 
the BDMPS-Z formalism, but takes into account finite- 
length effects. These effects include vacuuum-medium 
interference effects, the formation time effect at small 
X and the large-angle cut-off k < k±, although these as- 
pects are not easily separated in the calculation [53] . The 
effect of the finite-length corrections is a pronounced re- 
duction of the radiation at small gluon energies k, which 
is probably mostly due to the large angle cut-off. 

The green dashed curve shows the AMY result with 
the large-angle cut-off described in the previous section. 
This reduces the radiation at small k, as expected, but 
the reduction is not as pronounced as in the ASW-MS 
calculation (compared to BMDPS-Z). The smaller effect 
of the cut-off in AMY is likely due to the softer, Gaussian, 
cut-off that is used in AMY, compared to a hard cut-off 
k± < k in ASW-MS, although it is not excluded that 
the transverse momentum spectra are different, which 
would also affect the cut-off. The dip in the AMY result 
around k — E is caused by the cutoff (37) for the case 
that a gluon with momentum k ~ p is radiated. In that 
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FIG. 13: (Color online) Ratio of the rate for the process q — > qg suppressing emission angles larger than n/2 (left panel) and 
7r/4 (right panel) and the original AMY rate for two different a s and p = 40 T. The lines for g — » gg and q — > qg coincide. 




FIG. 14: (Color online) Gluon radiation spectra as a function of gluon momentum k for AMY (with and without fcx,max 
cut-off), ASW-MS, and the BDMPS-Z multiple soft scattering result Eq. 30 (with and without small x approximation) for a 
medium of length L = 2 fm (left panel) and L = 5 fm (right panel) . 



case the outgoing quark's longitudinal momentum p — k 
is very small and even a small kj> can lead to a large 
angle for the outgoing quark. Hence, the cutoff becomes 
relevant and suppresses the rate in this case as it does 
for small outgoing gluon momentum. 



D. The Higher-Twist (HT) approach 

In the application of perturbative QCD to hard pro- 
cesses, a fundamental role is played by the so called fac- 
torization theorems [ST]. These theorems demonstrate 
that all higher order corrections to a particular hard pro- 
cess may be cast in a form where the hard, short-distance 
part of a particular process can be separated order-by- 
order in coupling from the long-distance, soft parts of the 
process as long as one ignores terms which are suppressed 
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by powers of the hard scale of the process (generically de- 
noted as Q 2 )- In the case of the simplest hard process: 
The total cross section in the deep-inelastic scattering 
(DIS) of an electron on a proton, these power corrections 
can be expressed in the simpler language of the Opera- 
tor Product Expansion (OPE). In the OPE, the singular 
product of operators at closely separated points can be 
expressed as a series of terms which contain the product 
of progressively more complicated local operators and 
progressively less singular c-functions. The power sup- 
pression (i.e. power of Q 2 in the denominator) of each 
term depends on the "twist" of the local operator defined 
as the difference of the energy dimensions of the operator 
and the highest spin that can be constructed through the 
product of the various individual operators in the local 
product. 

So far, factorization has been rigorously proven for the 
leading twist part of a set of high Q 2 processes such as sin- 
gle hadron inclusive e + e~ annihilation, Drell-Yan, DIS, 
hadron-hadron collisions etc. Factorization at next-to- 
leading twist has been proved for a sub-class of processes, 
such as the total cross section in DIS [73]. As a result, 
there now exists a formulation of a class of higher twist 
corrections to the total cross section in DIS [75 . The 
higher twist energy loss approach is based on applying 
this formalism to the problem of calculating the change 
in the fragmentation function of a quark produced by the 
DIS of an electron on a nucleus [49] [76] . While the initial 
attempt focused on computing solely the single scatter- 
ing correction to the one gluon emission cross section, 
this has been subsequently generalized to include both 
multiple emissions in a DGLAP like formalism [50], as 
well as multiple scattering [77] . In the following we will 
describe the pertinent details of the HT approach assum- 
ing that a quark is produced in a hard collision at one 
edge of the QGP brick. We will first describe the setup 
where a quark is produced in the DIS of an electron on a 
nucleus and then having factorized the initial state and 
hard cross section from the final state evolution of the 
produced quark, we will consider the change in the evo- 
lution of this quark as it moves through the brick. 



1. Gluon emission formalism 

In this section, we describe the radiative part of the 
higher-twist (HT) calculation as applied to the "brick 
problem". In short, the higher- twist calculation consists 
of including a class of medium corrections to the process 
of jet evolution in vacuum, brought about by the mul- 
tiple scattering of the hard partons in a medium. It is 
most straightforwardly derived in the case of single inclu- 
sive deep-inelastic scattering in a large nucleus, with the 
nucleus playing the role of the medium. While most scat- 
tering corrections are suppressed by powers of the hard 
scale Q 2 , a subset of these are enhanced by the length 
of the medium and these are included in the calculation. 
Thus, the expansion parameter in the HT approach is 



a s qL/Q 2 where q is the transverse momentum squared 
imparted to a single parton per unit length and L is the 
length traversed by the parton. 

Consider the case of deep inelastic scattering (DIS) 
on a nucleus (in the Breit frame). The nucleus has 
a large momentum in the positive light cone direction 
0, 0, 0] with p + the mean momentum of a nucleon. 
The incoming virtual photon has a momentum which 
may in general be expressed as [— Q 2 /2q~, q~, 0, 0]; in 
the Breit frame q~ = Q/y/2 (light-cone coordinates in 
this section use the convention p + = ^(p + P 3 ) and 

p~ = -^{p° — P 3 ))- The inclusive cross section to pro- 
duce a hard hadron, which carries a momentum fraction 
z of the initial produced hard quark may be expressed in 
a factorized form as 



da 

dz 



dxG(x,Q 2 )~D(z,Q 2 ), (38) 



where G(x) is the parton distribution function (PDF) 
to obtain a hard quark in the nucleon with momentum 
fraction x. In the Breit frame the momentum of the in- 
coming quark is xp + — Q/y/2. Thus the produced quark 
has an outgoing momentum q~ = Q/y/2. The produced 
quark is virtual with a virtuality smaller than the hard 
scale usually denoted as XQ where A « 1. The other 
two factors are the hard partonic cross section da/dQ 2 
and the final fragmentation function D(z, Q 2 ). The scale 
in the fragmentation function is the factorization scale 
and also represents the maximum possible virtuality of 
the produced hard jet. The fragmentation function at 
the scale Q may be obtained from a measured fragmenta- 
tion function at a lower scale using the DGLAP evolution 
equations 



dD h q { Zl Q 2 ) = q s (Q2) fUy 
91og(Q 2 ) 2ir 



There is an implied sum over Savors i which includes 
all possible types of partons that may split off from the 
hard leading quark denoted as q. The kinematics of the 
scattering and emission process may be illustrated with 
the diagram in Fig. [15] 

In the case of DIS on a large nucleus, the above fac- 
torized form may be assumed to hold with the only 
change being the replacement of the vacuum evolved frag- 
mentation function with a medium modified fragmenta- 
tion function (as well as a replacement of the nucleon 
PDF with a nuclear PDF). The medium modified frag- 
mentation function contains the usual vacuum evolution 



piece Eq. ( 39 ) and a medium piece which includes both 



terms which are interferences between medium induced 
radiation and vacuum radiation as well as terms where 
both the amplitude and the complex conjugate repre- 
sent medium induced radiation. Once so factorized, the 
medium modified fragmentation function can be used to 
compute the single hadron inclusive cross section in any 
process by simply replacing the initial state parton dis- 
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4 

FIG. 15: A scattering and emission diagram with the cor- 
responding momenta, using the notation convention of the 
Higher Twist formalism. 



tribution and hard cross section by those appropriate for 
the process in question. 



For the Brick Challenge we ignore all the initial state 
functions and hard cross sections. We assume that the 
quark is produced at one edge of the brick designated 
as the origin and travels in the negative light cone direc- 
tion with a negative light cone momentum q~ . We assign 
the quark an initial virtuality Q 2 . Since this is not the 
Breit frame in DIS there is no implied relation between 
q~ and Q 2 . The equation for the medium modified frag- 
mentation function with an initial light-cone momentum 
q~ and virtuality Q 2 , which starts at the location £~~ and 
exits at £7, is given as 



dD%(z,Q 2 g-)\% = as 
dlog(Q 2 ) - 2tt 



i Cf 
dy 

y 

Ci 



dtP{y)K q -, Q2 {y,QD h q 



-,Q 2 ,q y 



y 



c 



(40) 



In the equation above, we have dropped the light-cone (— ) superscript on the positions. Note that the medium 
modified fragmentation function is now not only a function of Q 2 and z but also of q~ and £. The calculation of 
the evolution equation then requires the evolution of a three dimensional matrix (in The medium kernel 

Kq-,Q 2 (y, f° r a quark jet is given as |77j 



Kq-,Q'(y,t) = 



(q(C)-(l-y)l + (l-y) 2 q Q ) 



In the equation above, q(() without any subscripts rep- 
resents the position (£) dependent transport coefficient 
of a gluon, which may be expressed in operator form as 



8ir 2 a a CA 
N 2 - 1 



dy-(X\Tr [f^(»-KO<r] l*X42) 



where, F a ^ v is the position dependent field strength ten- 
sor of the gluon field and y~ represents the light-cone 
separation between these two field insertions. The state 
\X) represents the matter through which the jet propa- 
gates. We note that in the evaluation of multiple scat- 
tering diagrams leading to Eq. [41] it is assumed that 
the transverse momentum exchanged with the medium 
kj_ is soft, so that a Taylor expansion can be made. The 
leading non-zero term in this expansion is the term pro- 
portional to the square i.e. |fcj_| 2 - In combination with 
the gluon vector potential this term yields the transport 
coefficient q. The effect of higher terms in the Taylor 



2 - 2 cos 



QHt-d) 
29-1/(1 - y) 



(41) 



r 



expansion of k± is ignored, so that distribution of k± 
is approximated as a two-dimensional Gaussian with a 
mean k± = and a width proportional to q. 

The q for a quark scattering off the gluon field is related 
to the above expression as q~Q = In the factorized 

formalism of the HT approach, there is no scheme to eval- 
uate these coefficients and thus their values are external 
to the formalism. In most cases these are taken as fit 
parameters which are dialed to match one data point. 



2. Single gluon emission 



A gluon emission spectrum dN/dz can be calculated 



by integrating Eq. 40 over Q 2 using the gluon radiation 
kernel Eq. [41] and the gluon momentum fraction z 
giving 
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20 -r 



dN/dz 10 




on the radiated gluon momentum fraction, 



FIG. 16: The single gluon emission kernel in the HT approach 
for a quark with energy E — 20 GeV in a static medium of 
length 4 fm which has a <jq = 1 GeV 2 /fm with the maximum 
momentum cutoff from the medium set at M = 1.2 GeV 
(black solid line). The variable z here is equivalent to the 
variable x + used in the WHDG approach, earlier in this paper. 
Also shown are the single gluon spectra for M = 2.4 GeV (red 
dashed line) and M = 0.6 GeV (green dot-dashed line). See 
text for details. 



where we have taken only the first term of the emission 
kernel, and ignored terms which are suppressed by fac- 
tors z = (1 — y) (soft gluon approximation). The lower 
bound for the Q 2 integral is given by the requirement 
that the formation time of the gluon should be smaller 
than the length of the medium: Q^in = E/L. The upper 
integration limit is given by 



Ql 



2/(1 - y)2q-p+ = z(l - z)2ME, 



(44) 



which is discussed further below (cf. Eq. 45 1 



The medium induced single gluon emission distribution 
in momentum fraction z is plotted in Fig. |16[ for a 20 
GeV quark, propagating through a 4 f m long static brick, 
using M = 3T. 

The medium induced gluon emission kernel has soft di- 
vergences at y = 0, 1 as is the case for the gluon emission 
kernel in vacuum. However, given that q is very similar 
to a gluon distribution function it depends on the mo- 
mentum fraction of the gluon off which the hard parton 
scatters. In the case of DIS, this momentum fraction, 
referred to as xl, is given as xl = Q 2 /{2q~p + y(l — y))- 
Where, p + = M/y/2, M is the mass of the nucleon and 
q~ = \/2E, where E is the energy of the virtual pho- 
ton in the nucleus rest frame. One expects q{xi,) to be 
dominated by small values of the fraction xl and rapidly 
vanish for values of Xl > 1. This introduces a constraint 



< V < 




Q 2 



2q-p^ 



(45) 



Note that in the above equation Q 2 ~ k\ the transverse 
momentum of the radiated gluon. Since only values of 
y which lie between and 1 are allowed, the value of 
k\ is always constrained to be less than q~p + where p + 
is the momentum of the proton. In a multiple scatter- 
ing scenario, this constraint only applies to the largest 
k±, i.e. the first emission. Successive emissions are or- 
dered in k±. While, this limit is well defined in DIS, it 
is slightly uncertain in the case of a quark gluon plasma 
created in a heavy-ion collision. However, one argues 
that the probability of a medium at a temperature T to 
radiate a gluon with momentum much larger than the 
mean momentum of 3T should be very small. Thus one 
replaces p + = M/V2 — > 3T/^/2 in the equation above. 
Here, M represents the largest momentum that may be 
exchanged between the medium and the jet, which is 
achieved when a thermal parton from the medium di- 
rectly impacts with the jet. The central black solid line 
in Fig. [l6|corresponds to the choice of T = 400 MeV, with 
M = 3T = 1.2 GeV. The red dashed line corresponds to 
the replacing the upper limit with M = 6T = 2.4 GeV 
and the dot-dashed green line to replacing the upper limit 
with M = 3T/2 = 0.6 GeV. We note that the main dif- 
ference occurs at at low z, i.e. in the soft gluon limit, as 
would be expected 



3. Multiple gluon emission 

We now plot the result of multiple emissions from the 
single quark parton produced at one edge of a static brick. 
In the soft gluon limit, the medium modified fragmenta- 
tion function of a quark is obtained by replacing the vac- 
uum splitting function with the sum of the vacuum split- 
ting function and the medium modified splitting function 
in Eq. ( 40 1 . These equations require an input fragmen- 



tation function at some lower initial scale usually chosen 
to be 1 GeV. The input fragmentation function is usually 
chosen to be the vacuum fragmentation function. Note 
that even if the initial input fragmentation function is 
assumed to have no position dependence, the evolution 
equation generates a position dependence in the evolved 
fragmentation functions at any higher scale. A more cor- 
rect input is to have a position dependent fragmentation 
function which is the vacuum fragmentation function at 
distances beyond the brick and goes to zero swiftly for 
fragmentation within the brick. For the HT plots in this 
paper we will ignore this sophistication and ignore the po- 
sition dependence of the fragmentation function. Even at 
this level of approximation, the equations above are far 
too numerically intensive to solve. One usually replaces 
the position dependence with the initial position C, — > Q. 
The evolution equations now represent the evolution of 
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q = 1 GeV 2 /fm, E = 20 GeV, L = 4 fm 



l ■_ 




FIG. 17: The results of evolving an input [q(q) — > q(q)] frag- 
mentation function at fj 2 = 1 GeV 2 given by a 6(1 



(see 

Eq. ([46}) up to a higher scale of Q 2 = 10 GeV 2 in a 4 fm 
brick. The transport coefficient <jq = 1 GeV 2 /fm with the 
maximum momentum cutoff from the medium set at 1 GeV 
(solid black line), 2 GeV (red dashed line) and 0.5 GeV (blue 
dash-dotted line). See text for details. 



a two dimensional matrix and these represent the calcu- 
lations which will be presented in this paper. A further 
approximation which is sometimes used in the literature 
is to also drop the energy dependence of the fragmenta- 
tion functions. This approximation will not be made and 
the energy dependence of the medium modified fragmen- 
tation function will be retained. 

In order to compare with other formalisms which 
present the energy loss of a single quark, we introduce 
a quark-to-quark fragmentation function as the input at 
the lower scale of /i 2 = 1 GeV 2 . In actuality our input 
is a singlet quark fragmentation function ( (q + q)/2 — > 
(q + q)/2). Our input at 1 GeV may be expressed as 

«±a(z, lGeV 2 ) = 8(1 -z) 
LW^„±=(z,lGeV 2 ) = 
D^s+aiz, lGeV 2 ) = 0. 



(46) 



In Fig. [T7J these are represented by the solid green line 
[u(uj — ¥ u(u)] along the z = 1 axis and the cyan dot- 
dashed line [g — > u(u)] along the dN/dz = axis. We 
do not plot the results of the d(d) —¥ u(u) fragmentation. 
These distributions are then evolved up to the higher 
scale of Q 2 = 10 GeV 2 . The results for the three different 
momentum bounds from the medium M — 1.5T, 3T and 
6T (as discussed above) with T = 1[3 GeV and a quark 
4q = 1 GeV 2 /fm are plotted in Fig. 
quark and gluon distributions are 
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Note that the full 
evolved. This means 
that if in one step of the evolution a hard quark loses 



energy by splitting into a quark and a gluon, in the next 
step of the evolution both the softer quark and gluon will 
further lose energy by splitting themselves. 

As would be expected, we note that in the case when 
the bound from the medium is set at the higher value 
of M — 2 GeV, there is more gluon radiation and thus 
there is more modification of the input distribution and 
vice-versa for the case when the bound is set at a lower 
value of M = 0.5 GeV. The dependence of the evolved 
distribution on the bound is enhanced primarily due to 
the numerical instability of the vacuum DGLAP evolu- 
tion equations on an input distribution that is singular at 
z — > 1. Note that the evolution equations include a pure 
vacuum part and a medium induced part as expressed 
in Eq. (40). The pure vacuum splitting functions are 



defined using the (+)-functions, i.e., 

dD(z,Q 2 ) ^ f 1 dy (z 



(47) 



- f ld ^P(y) D ( Z -,Q 2 )^D(z,Q 2 ) fdyP(y). 

n Jz v \y J Jo 



The above equation assumes that D(z) is well defined at 
all values of z over which the evolution is carried out. 
This condition is not fulfilled by the strict <5-function. 
The calculations presented in Fig. [17] are carried out by 
using a regulator to define the <5-function. We should 
point out that only one choice of regulator has yielded 
convergent results as the regulator is removed from the 
results, viz. 



5(1 -z) = 



1 



V z > 1- e, 



0, V z < 1- e, 



(48) 



in the limit e — > 0. Other choices of regulation, e.g. a 
narrow Gaussian with a width given by e turned out to 
yield results which were non-convergent within the limits 
of numerical accuracy used. 



4- Medium modified fragmentation function 

The DGLAP evolution equations used by the HT 
scheme require that the input fragmentation function be 
well defined at all values of momentum fraction z. To 
illustrate the degree of control in these calculations we 
now repeat the calculation above but using a well de- 
fined input distribution: a tt° fragmentation function (in 
this case taken from the KKP parametrization [75]). The 
input is taken at /j, 2 = 1 GeV 2 and evolved to 10 GeV 2 . 
The results are plotted in Fig. [18] The brick is again of 
length 4 fm and has a quark <?q = 1 GeV 2 /fm. Since 
fragmentation functions vary over orders of magnitude 
as z runs from to 1 we plot the ratio of the medium 
modified fragmentation function to a vacuum fragmen- 
tation function at the same scale. The denominator is 
obtained by evolving the fragmentation function from 1 
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FIG. 18: The ratio of the medium modified fragmentation 
function (of a ir° from a quark with E = 20 GeV) to the 
vacuum fragmentation function at Q 2 — 10 GeV 2 . The 
input in both cases is the KKP fragmentation function at 
fi 2 — 1 GeV 2 . The medium is a static brick with a length of 
4 fm and a qQ = 1 GeV 2 /fm. Three different values for the 
maximum momentum cutoff from the medium are used 1 GeV 
(black solid line), 2 GeV (red dashed line), and 0.5 GeV (blue 
dot-dashed line). The open red squares and open blue circles 
indicate the results with qg changed to make the lines with 
M — 2 GeV and M = 0.5 GeV approximately coincide with 
the central line with M — 1 GeV. See text for details. 



GeV 2 to 10 GeV 2 with a q Q = 0. The black solid line 
represents the standard case with a bound on the mo- 
mentum from the medium set at M — lGeV. The red 
dashed and the blue dot-dashed lines represent the cases 
with M = 2 GeV and M — 0.5 GeV. As expected, raising 
the bound, increases the phase space of gluon radiation 
and thus increases the modification. 



5. Uncertainty due to large angle radiation 

The uncertainty in the bound on the momentum that 
may emanate from the medium introduces an uncertainty 
in the extracted value of q. To estimate this uncertainty, 
we take the case with M = 0.5 GeV (the blue dot-dashed) 
which shows less modification than the case with M = 1 
GeV and increase the q until the blue dot-dashed line 
coincides with the black line. The final result is indi- 
cated with a blue dot-dashed line with open circles. This 
requires a 30% increase in qq. Similarly we reduce the 
qQ for the case with M = 2 GeV (red dashed line) to 
make it coincide with the black line. The final result is 
indicated by the red dashed line with open squares and 
requires a qQ which is 10% lower than the central value of 
1 GeV 2 /fm. Thus the uncertainty in the extracted value 



of q in HT calculations is between 10% to 30% when the 
qQ ~ 1 GeV 2 /fm and the length of the medium is of the 
order of 4 fm. The reader will note that these values are 
quite representative of the respective values at RHIC and 
thus provide the range of error in estimates of q extracted 
in RHIC collisions. 



IV. SYSTEMATIC COMPARISONS 

In this section, the multiple-soft scattering approxima- 
tion is compared with the opacity expansion formalism 
and the McGill AMY implementation. Specifically, the 
following energy loss models are compared: 

• ASW-MS: The multiple-soft scattering approxima- 
tion as formulated by ASW [33] and discussed in 
Sections HIA] and HICl 

• ASW-SH: The single hard scattering approxima- 
tion as described formulated by ASW [39] . and 
shown in Eq. ([3] and Table [I] The original formu- 
lation uses a fixed value for L/X = 1 for analytical 
convenience. For this work L/X is calculated from 
the temperature T in the medium (see Appendix 

• DGLV: The single hard scattering approximation 
using the same choices for in-medium gluon and 
quark masses and the kinematic bounds as used by 
WHDG @7] , as given in Eq. |] and Table [I] 

• AMY radiative: The AMY formalism based on a 
thermal effective field theory as described in Sec- 
tion |III C| In the interest of the comparison to 
the other formalisms, the gluon splitting process 
9 99 is removed from the calculation. 



A. Comparison at fixed temperature 

Figure [20] shows the single gluon energy spectra for 
medium-induced radiation for a uniform medium with 
T = 300 MeV and two different path lengths L = 2 fm 
(left panel) and L = 5 fm (right panel) using the four dif- 
ferent formalisms discussed above. The HTL-based AMY 
formalism (green curve) used here is without the large- 
angle cut-off discussed in Section III C 3 and it gives the 
largest amount of radiation. The opacity expansions (red 
curves, DGLV and ASW-SH) and the Higher Twist for- 
malism (black curves, HT) are next while the multiple- 
soft scattering approximation ASW-MS gives the least 
radiation. For ASW-MS and HT two curves are shown, 
using different relations between the transport coefficient 
q and the temperature T. One is based on the soft- 
scattering limit from HTL field theory, Eq. 16 which 
gives q = 2CUa s T^ 2 log(g 2 nax //x 2 )) while the other is 
Eq. |A6[ based on Eq. |A5| The single gluon spectrum 
for the HT formalism is evaluated using Eq. 43 which 
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FIG. 19: (Color online) The single gluon distribution as function of gluon energy uj for a uniform medium with T = 300 MeV 
and two different lengths, L = 2 fm (left panel) and L = 5 fm (right panel). For the AMY calculation, the outgoing gluon 
spectra are plotted, including the evolution via the rate equation Eq. |21| The bands for the ASW-MS and HT formalisms 
indicate the uncertainty from evaluating q using Eq. |16| and |A6| 




FIG. 20: (Color online) The final quark energy spectrum as function oi xe = E out /E = 1 — e for a uniform medium with 
temperature T = 300 MeV and two path lengths L = 2 fm (left panel) and L = 5 fm (right panel). The symbols at xe — 
indicate the probability that a quark is absorbed and at xe = 1 the probability that a quark does not interact with the medium. 
Blue circles: ASW-MS. Open red squares: DGLV. red circles: ASW-SH. 



does not include modified DGLAP-evolution of the par- 
tons. The HT formalism gives more radiation at larger 
gluon energies than the opacity expansion and ASW-MS, 
while the peak at small oj is less pro nou nced for the de- 
fault Ql iax = 6z(l 
of changing 



z)ET. See Fig. 



16 



to see the effect 
t . We also note that the L = 5 fm case is 
beyond the strict validity of the HT formalism, since the 
lower virtuality bound Q^m = E/L sa 0.4 GeV 2 . Despite 



this, the result is shown for comparison. 

A similar ordering between the different formalisms is 
visible for L = 2 fm and L = 5 fm, although the differ- 
ence between the AMY result and the opacity expansions 
is larger at small L. The peak at small ui for the HT for- 
malism is also more pronounced for larger path length. 
The path-length dependence has been studied in more 
detail in [SHI, where is is shown that for small L, AMY 
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overestimates the radiation by ignoring finite-length ef- 
fects, while for large path-lengths, the first order (N = 1) 
opacity expansion overestimates the radiation by neglect- 
ing interference between subsequent scatterings. At large 
path-lengths, the AMY result is close to the ASW-MS re- 
sult, but only for large gluon energies. At smaller gluon 
energies, the large-angle cut-off strongly reduces radia- 
tion in ASW-MS. This effect is not implemented in AMY 



(see Section III C 3 1 



For ASW-MS and ASW-SH the thermal quark and 
gluon masses are set to zero while in DGLV those masses 
are finite. Including the thermal masses reduces soft 
gluon radiation which can be seen comparing the ASW- 
SH and DGLV curves at small w. Another difference 
as discussed in Section |IIIB| is the cut-off fc max on the 
transverse momentum of the radiated gluon. For ASW- 
MS and ASW-SH it is possible that a radiated gluon has 
larger momentum than the initial energy of the incom- 
ing quark because there is no kinematical bound on the 
gluon energy uj. 

To calculate the total energy loss, the gluon spectra 



from Fig. 19 are used as input for a calculation of mul- 
tiple gluon emission. For the opacity expansions and the 
ASW-MS calculation, a Poisson ansatz is used (Eq. [IJ, 
while AMY use rate equations (Eq 21 1. The resulting 

No re- 
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outgoing quark spectra are shown in Fig. 
suits for the HT formalism are shown, because medium- 
induced and vacuum radiation cannot be separated in 
the medium-modified DGLAP evolution, which makes it 
difficult to compare to the other formalism on the same 
footing. The Poisson ansatz used for the opacity expan- 
sions and ASW-MS gives a finite probability to the par- 
ton to lose no energy (zero gluons radiated) which is in- 
dicated by the symbols at xb = 1 in Fig. [20j At the 
same medium density and L = 5 fm, this probability is 
approximate twice as large for the multiple soft scatter- 
ing approximation as for the opacity expansions (right 
panel of Fig. 20). On the other hand, the gluons that 



are radiated are softer for the opacity expansions than 
for ASW-MS. For a brick of L = 2 fm and a tempera- 
ture of T = 300 MeV the discrete weights at %e = 1 are 
between 0.5 and 1 (see also Table III), so that the be- 



haviour is dominated by the discrete weights, although 
there are large probabilities for small energy loss xe ^ 0.8 
as well. The distribution from the AMY formalism is cut 
off at xe < 0.1, because the formalism is not valid when 
then parton energies approach the thermal energy; par- 
tons with energy less than 2 GeV are removed from the 
evolution. 



It is clear from Fig. 20 that the different energy loss 
formalisms do not result in similar outgoing quark distri- 
butions. 



B. Suppression factor in a QGP Brick 

To characterise the energy loss distributions in Fig. [20] 
in a single number, we calculate an approximation of the 



nuclear suppression factor Raa in the following way. 

The measured hadron spectra at RHIC approximately 
follow a power law: dN/dp? = Ap^ 11 ■ If the energy 
of each hadron is reduced by a fraction e, the hadron 
momentum after energy loss p' T = (1 — e)pr and the 
spectrum after energy loss 



dN 

dp'rp 



= A 



(l-e) n dpT =A (l-e)' 
pVj, dp' T p 1 ^ 



(49) 



In this case Raa = (1 — e)" _1 - This can be generalized 
to a distribution of energy loss, in which case the nu- 
clear modification factor Raa can be approximated by 
the appropriately weighted average over the energy loss 
probability distribution P(e) 



Rn — 



= [ de(l-ey 
Jo 



(50) 



with e = AE/E. In the following, Ri will be used as a 
proxy for Raa, because for RHIC energies the hadron px 
spectrum can be approximated by a power law spectrum 
with n = 6.5 for p T > 2. GeV/c (75]. 

Figure 21 shows the dependence of the suppression 
factor i?7 on q for the different formalisms. The fig- 
ure clearly shows that both opacity expansion formalisms 
generate a larger suppression at the same medium density 
(same q) than the multiple soft scattering approximation. 
The AMY formalism generates the largest suppression 
at a given density. The values of the transport coeffi- 
cient needed to reach a similar suppression as measured 
at RHIC HDHE3, R7 ~ 0.25, are listed in Table |H] and 
differ by a factor 5-10 between AMY and ASW-MS. 



C. Comparison at fixed suppression R7 

In Fig. [22] the inclusive gluon spectra are shown for a 
uniform medium of L — 2 and 5 fm, for a fixed suppres- 
sion i?7 = 0.25 using the medium density values given 
in Table |ll] The single gluon spectrum from DGLV does 
not extend beyond ui = E because that formalism im- 
plements a large momentum cut-off to impose forward 
propagation of the final state quark (cf. Eq. 11). For 
the AMY gluon spectrum only q — > q + g splittings are 



included. In the AMY formalism it is not possible to dis- 
tinguish between thermal and radiated gluons for u) < 2 
GeV which is why in this region for AMY the gluon spec- 
trum is not shown. We note that the ASW-MS single 
gluon spectrum at fixed suppression is harder than that 
obtained in the opacity expansions. 



Figure 23 shows the outgoing quark energy spectrum 
as function of xe — 1 — e for the two bricks of different 
lengths and -R7 = 0.25. The probability that a parton 
is absorbed in the medium is indicated by the markers 
at xe = 0. In this case the the energies of the multiple 
radiated gluons add up to a total energy loss that exceeds 
the initial energy of the parent parton and the parton is 
absorbed in the medium. The large x cut-off on the single 
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FIG. 21: Correlation between R7 and q for a primary quark with E = 20 GeV for different energy loss formalisms. The 
horizontal black dashed line indicates R7 — 0.25. 
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= 0.25 




T (MeV) 


q (GeV 2 /fm) 


ui c or uj c L/\ (GeV) 


7? or RL/X 


L/X 


m D (GeV) 






ASW-MS 


1030 


23.2 


236 


2393 






L = 


2 fm 


WHDG 


936 


17.8 


105 


1063 


6.25 


1.82 






ASW-SH 


727 


8.86 


49.0 


500 


4.85 


1.41 






AMY 


480 


2.7 














ASW-MS 


434 


2.11 


134 


3401 






L = 


5 fm 


WHDG 


358 


1.23 


36.5 


925 


5.97 


0.69 






ASW-SH 


326 


0.95 


27.6 


702 


5.44 


0.63 






AMY 


235 


0.4 











TABLE II: Values of the model parameters required to reach the typical suppression of Rj = 0.25. 



gluon spectra in the DGLV formalisms leads to a smaller 
probability for absorption than in ASW-MS and ASW- 
SH. The corresponding probability for losing no energy 
at all is given by the markers at Xe = 1- 

The probabilities for absorption and no interaction 
of the multiple-soft scattering approximation are larger 
than for the opacity expansions. The continuous part 
of the energy loss probability distribution is more rele- 
vant in the opacity expansion than for the multiple soft 
scattering approximation. 

In the AMY formalism, the outgoing quark spectrum 
peaks at xe = 0.9, which corresponds to an energy loss 
AE ~ 2 GeV. This radiation falls in the the region where 
we cannot distinguish the radiated gluons from the ther- 
mal gluons in the AMY formalism, as indicated by the 
cut in Fig. [22} 



1. Some characteristic quantities 

When all the models are tuned to a fixed amount of 
suppression i?7, the radiated gluon spectra are quite dif- 
ferent. To quantify the difference, we calculate a few 
characteristic quantities for energy loss which are tabu- 
lated in Table [TTT1 



The average number of emitted gluons can be obtained 
by integrating the gluon spectrum (e.g. Fig. 22 1 over the 
gluon energy u: 



(N g ) = / du 



dl 

dui 



(51) 



Note that this determines the probability for no energy 
loss: 



Po 



(52) 



The probability for complete absorption of the parton is 
given by 



i' 00 

Pl = / d(AE) P(AE), 
Je 



(53) 



where E is the energy with which the parent parton en- 
ters the brick. The total average energy loss for the in- 
coming parton is calculated from the energy loss proba- 
bility distribution P(AE): 



(AE) 



[ 'd(AE) P(AE), 
Jo 



(54) 



for which only the surviving partons are taken into ac- 
count. For the AMY formalism the integral starts at 
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FIG. 22: Inclusive gluon radiation spectrum for quarks with E — 20 GeV with in-medium path-lengths L = 2 fm (left panel) 
and L = 5 fm (right panel). The medium density has been tuned for each formalism to give i?7 = 0.25. The AMY result 
includes the evolution via the rate equations Eq. |21| 




FIG. 23: The final quark energy spectra as function of ie = 1 — e and for path-lengths L — 2 fm (left panel) and L — 5 
fm (right panel). All models are scaled to the same suppression Rr = 0.25. The squares at ie = indicate the probability 
that a quark is absorbed and at xe = 1 the probability that a quark does not interact with the medium. Solid blue squares: 
BDMPS-Z/ASW-MS. Open red squares: DGLV. Solid red squares: ASW-SH. 



negative energies (in principle — oo) because the parton 
can also absorb energy from the medium (see Figs 20 and 
231. 



CONCLUSION 



The values collected in Table [TIT] show that tuning the 
models to the same suppression factor Rj does not imply 
that the mean energy loss is the same. On the other hand, 
when the suppression factor Rj is fixed, the number of 
radiated gluons is similar. 



For a summary of the main findings we refer back to 
Section [ETDj 
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Po 


Pi 






AMY 




3.95 








6.58 








— 2 f m 


ASW-MS 


0.17 


0.3 


0.84 


0.01 


1.58 


5.9 


0.2 


0.46 






WHDG 


0.48 


1.24 


0.62 


0.00 


1.94 


7.7 


0.14 


0.14 






ASW-SH 


0.25 


0.93 


0.78 


0.03 


1.92 


6.9 


0.15 


0.35 






AMY 




8.30 








6.47 






L 


= 5 fm 


ASW-MS 


1.17 


3.95 


0.31 


0.08 


1.83 


6.4 


0.16 


0.33 






WHDG 


1.96 


5.63 


0.14 


0.05 


2.44 


7.3 


0.09 


0.11 






ASW-SH 


2.09 


5.76 


0.12 


0.2 


2.47 


6.7 


0.08 


0.26 



TABLE III: Table of different variables for the different energy loss models. Left columns: for T = 300 MeV for all models. 
Right columns: corresponding medium properties required to reach the typical suppression Rr = 0.25. 
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Appendix A: Common definition of medium 
properties q, fi and A 

In this Appendix we review a number of commonly 
used equations to relate properties of the medium, such 
as the transport coeffient q, mean free path A and tem- 
perature T. 



1. The basis: a HTL plasma 

The basis for all the equations used are free scatter- 
ing cross sections, supplemented with elements of Hard 
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Thermal Loop (HTL) field theory to screen the soft di- 
vergence in the free cross section. 

For a medium with temperature T, the Debye screen- 
ing mass mo is 



l 2 D = (1 + -N f )g 2 T 2 



(Al) 



where g is the coupling constant (g = 47ra5) and Nf is 
the number of quark flavours (0 for a pure gluon gas). 
We also define a number density of the medium Af |83j : 



^l(l+ 1 N f )T 3 . 
C(2) V 4 /; 



(A2) 



The differential scattering rate for a hard parton 
traversing the medium is known for two limits of the 
exchanged momentum q± |83j : 



9 2 Tm 2 D 



d 2 q±_ (2tt) 



qi(ql 

IK 



if |qi| «T 



if |qj »T, 



(A3) 



the relevant color factor Cr for most energy loss calcula- 
tions is the color factor of the gluon Cr = Ca = 3. 



2. Transport coefficient q 

The transport coefficient q, which is the only medium 
parameter that enters in the multiple soft scattering ap- 
proximation, is defined as the mean momentum kick 
squared per unit path length 



P 



i2 2 

a qrq T 



da 
d 2 q T 



,2 dT e i 2 
d q±-^ — q ± , 
d 2 q±_ 



(A4) 



in which T e i is the rate for elastic collisions in the plasma, 
qj^ is the transverse momentum exchanged with the 
medium g max is the upper limit for the exchanged mo- 
mentum, which represents an ultraviolet cut-off. 

The transport coefficient q integrates over all momen- 
tum transfers q and can therefore in general not be based 
on only the low or the high q± limit given in Eq. |A3| The 
following expression provides a smooth interpolation be- 
tween the two limits 



dT„ 



Cr 



5 4 A/- 



d*q x ~ (2tt) 2 qi(qi 



(A5) 



This interpolation is il lustr ated in Fig. [24] which shows 
low energy limit of Eq. ( A3 ) (curve labeled Soft) together 



with Eq. A5 (curve labeled Interpolated). In the soft re- 
gion the interpolated curve deviates by 17% due to the 
difference in the numerator and approaches the high en- 
ergy limit smoothly at large q±. 

The resulting transport coefficient is 



q(T) 



C R g A N(T) 

47T 



hi 



D 



(T) 



i 



(A6) 
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FIG. 24: Differential elastic cross-section as function of q± 
for the different energy limits. The curves labeled with In- 
terpolated correspond to the expression given in Eq. \A5\ . 
Curves labeled with Soft correspond to the low energy limit 
in Eq. (|A3|). 
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FIG. 25: The transport coefficient as function of the energy 
of the parton traversing a medium of different tempertures. 
The transport coefficient is calculated in the hard thermal 
loop formalism as in equation | A6| 



In the high energy limit the constant in the argument 
of the logarithm vanishes but in the low energy limit its 
presence is crucial because otherwise unphysical negative 
values for q are possible. In a realistic medium created in 
a heavy-ion collision naturally both cases will be present. 

The value for g max is expected to fall in the hard scat- 
tering regime, where the scattering centers can be ap- 
proximated as static, leading to <7 max = g(ET 3 ) 1 ^ 4 as 
has been argued in [73] • The dependence of q max on 
the energy of the incoming parton E leads to a loga- 
rithmic dependence of q on E. Figure [25] compares the 
resulting energy-dependent q to a commonly used energy- 
independent expression for q 



(A7) 
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where it is assumed that the mean squared momentum 
exchange per scattering is equal to mjj and the mean free 
path A is calculated as given below. 



3. Mean free path 



which is very similar to Eq. ~KE except for a change in 
the denominator which is needed to obtain a finite result 
for A. 

This leads to the following expression for the number 
of scattering centers 



The mean free path A used in the opacity expansion is 
usually calculated based on a Quark Gluon Plasma with 
static scattering centers, using the scattering rate 



Cr 
(2tt)2 



5 W 



(qi 



^2 \2 ' 
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